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Uncertainty  and  time-delay  in  real  systems  constitute  the  two  major  challenges 
that  face  control  engineers  since  both  can  contribute  to  instability  or  poor  perfor- 
mance. In  this  dissertation  three  analysis  and  control  design  problems  are  addressed. 
These  problems  involve  linear  systems  with  parametric  uncertainty  structure,  and 
linear  and  bilinear  systems  with  time  delay. 

In  the  first  problem,  the  Nyquist  robust  sensitivity  margin  is  proposed  as  a 
scalar  metric  for  robust  stability  and  robust  performance.  The  work  was  motivated 
by  the  critical  direction  theory  (CDT)  in  which  attention  was  given  to  plants  that 
lie  along  the  critical  direction.  The  advantage  of  the  new  metric,  however,  is  that  it 
takes  into  account  plants  that  are  in  close  proximity  to  the  critical  point  — 1 -h  jO  but 
that  do  not  lie  along  the  critical  direction.  The  approach  introduced  has  therefore 
the  advantage  of  capturing  the  worst  case  sensitivity  as  well  as  providing  a more 
meaningful  indication  of  robust  stability.  The  concept  has  been  applied  successfully 
to  a class  of  linear  systems  with  affine  uncertainty  structure. 


IX 


The  second  problem  involves  designing  a sliding  mode  control  (SMC)  to  stabilize 
a class  of  time-delay  linear  systems.  The  delay  is  assumed  to  exist  in  both  the  control 
variable  as  well  as  in  the  state  vector.  The  system  is  first  rendered  input-delay  free 
through  an  appropriate  transformation.  Then  an  SMC  is  designed  for  the  state- 
delay  system.  Sufficient  stability  conditions  ensuring  the  asymptotic  stability  of  the 
closed-loop  system  have  been  derived. 

The  third  problem  addresses  the  stabilization  of  a class  of  time-delay  bilinear 
systems.  A state- feedback  control  law  is  designed  to  ensure  the  asymptotic  stability 
of  the  delayed  bilinear  system.  The  work  builds  on  two  simple  scalar  systems  and 
utilizes  the  results  to  prove  a more  complicated  system.  The  analysis  allowed  us  to 
obtain  a bound  on  the  maximum  value  of  the  delay  that  the  system  can  tolerate. 
Furthermore,  a region  of  attraction  based  on  the  initial  condition  of  the  systems 
states  is  established. 
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CHAPTER  1 
INTRODUCTION 


1.1  Robustness  Analysis 

The  robustness  analysis  problem  investigates  the  behavior  of  a dynamical  system 
under  uncertainty,  namely,  how  the  system  stability  and  performance  are  influenced 
by  the  uncertainty.  Many  robust  stability  tools  have  been  developed  along  the  years, 
among  which  are  the  well-known  scalar  stability  margins:  the  structured  singular 
value  introduced  by  Doyle  [17]  and  the  multivariable  stability  margin  given  by 
Safanov  km{u))  [46]. 

The  critical  direction  theory  introduced  by  Latchman  and  Crisalle  [35]  and  later 
generalized  by  Baab  et  al.  [2]  also  provides  an  effective  tool  for  analyzing  the  robust 
stability  of  uncertain  systems,  namely,  the  Nyquist  robust  stability  margin  ki^{uj). 
The  concept  was  applied  successfully  to  a class  of  linear  systems  with  affine  and 
ellipsoidal  uncertainty  structure,  and  it  works  for  the  case  of  convex  and  non-convex 
value  sets. 

Uncertainties  are  classified  depending  on  their  source  as  non-parametric  (un- 
structured) and  parametric  (structured)  [7].  The  non-parametric  uncertainties  do 
not  have  a well  defined  structure  and  are  represented  by  a disk  which  over  bounds 
the  actual  uncertainty.  Therefore,  this  type  of  uncertainty  description  usually  intro- 
duces conservatism.  Examples  of  uncertainties  that  are  represented  as  unstructured 
include  non-linearities  and  un-modelled  dynamics.  Parametric  uncertainties,  on  the 
other  hand,  have  a structure  that  reflects  the  variation  of  the  system  parameters. 
Thus,  they  are  less  conservative.  Examples  of  such  uncertainties  include  interval  and 
ellipsoidal  uncertainty. 
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1.2  Sliding  Mode  Control 

A variable  structure  system  (VSS)  is  a dynamical  system  composed  of  distinct 
structures.  A VSS  switches  between  the  different  structures  based  on  the  value  of 
its  states  and  according  to  a switching  logic  which  takes  into  account  the  desired 
properties  in  each  structure.  In  fact,  a variable  structure  system  can  have  properties 
that  are  not  existent  in  its  individual  structures  [51]. 

A sliding  mode  control  system  (SMC)  is  a specific  case  of  VSS  in  which  the 
system  trajectories  exhibit  a sliding  behavior.  The  design  of  an  SMC  consists  of 
two  stages.  The  first  stage  is  the  design  of  a switching  surface  such  that  once  the 
trajectories  are  confined  to  the  surface  the  system  demonstrates  the  desired  properties 
(be.,  tracking,  regulation,  etc.)  The  second  stage  involves  the  design  of  a control  law 
that  forces  the  trajectories  into  the  sliding  manifold  (discontinuous  control),  and  a 
linear  feedback  control  that  guarantees  closed-loop  stability  (equivalent  control).  The 
latter  is  derived  by  setting  the  time  derivative  of  the  switching  function  equal  to  zero 
and  solving  for  the  control  law.  The  former  is  proposed  with  appropriate  gains  to 
allow  the  system  to  overcome  uncertainties. 

The  system  motion  in  SMC  runs  through  two  phases.  The  first  phase  (reaching 
phase)  is  characterized  by  a fast  motion.  The  system  during  this  phase  is  robust 
against  uncertainties  (matched  and  unmatched)  and  external  disturbances.  This 
is  mainly  due  to  the  discontinuous  control  law  which  acts  as  a high  gain  feedback 
control  that  counteracts  high  frequency  signals.  The  second  phase  (sliding  phase)  is 
characterized  by  a slow  motion.  The  system  is  however  robust  only  against  matched 
uncertainty. 

The  theory  of  sliding  mode  control  has  been  covered  comprehensively  in  the 
literature.  Utkin  [51]  presents  a survey  for  the  early  contributions  in  SMC.  The 
survey  by  Hung  et  al.  [26]  presents  a tutorial-like  paper  for  variable  structure  control 
(VSC)  with  sliding  mode.  An  interesting  tutorial  paper  by  DeCarlo  et  al.  [16] 
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provides  an  introduction  to  variable  structure  control  for  multivariable  nonlinear 
time-varying  systems.  Finally,  a useful  guide  to  SMC  is  also  given  by  Young  et  al. 
[56]. 

1.3  Bilinear  Systems 

Bilinear  systems  occupy  an  intermediate  level  between  linear  and  nonlinear  sys- 
tems in  terms  of  their  complexity,  the  general  form  of  a bilinear  systems  is  given 
as 

x{t)  = Ax{t)  + Bu{t)  + Nx{t)u{t)  (1.1) 

where  it  is  clear  that  the  control  action  enters  the  system  linearly  through  the  term 
{Bu{t)),  and  nonlinearly  through  the  term  (Nx{t)u{t)),  hence  the  name  bilinear 
system.  A special  form  of  the  system  (1.1)  is  the  homogenous  bilinear  system 

x{t)  = Ax{t)  A Nx{t)u{t) 

where  the  linear  part  is  omitted.  A formal  definition  of  a bilinear  system  is  given 
in  Elliott  [19].  Many  natural  as  well  as  man-made  systems  can  be  represented  as 
bilinear  models  [19,  39,  40]  . Examples  of  bilinear  systems  can  be  found  in  economics, 
industrial  processes,  and  biochemistry,  just  to  mention  a few. 

1.4  Thesis  Structure 

The  thesis  is  organized  as  follows.  In  Chapter  2,  a new  metric  to  the  robust 
stability  of  closed-loop  systems  with  affine  uncertainty  structure  is  presented.  The 
new  concept  is  motivated  by  the  fact  that  the  critical  direction  theory  considers  only 
plants  that  lie  along  the  critical  direction  defined  as  the  ray  starting  at  the  nominal 
plant  and  pointing  towards  the  critical  point  —1  + jO.  Hence,  plants  that  are  very 
close  to  the  critical  point  but  that  do  not  lie  on  the  critical  ray  are  ignored.  Therefore, 
the  Nyquist  robust  sensitivity  margin,  k^^s,  is  proposed  to  take  into  accounts  such 
plants.  Chapter  3 considers  the  stabilization  of  a class  of  time-delay  linear  systems 
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via  sliding  mode  control.  The  delayed  system  is  assumed  to  have  a constant  delay 
in  both  the  input  and  the  state.  In  Chapter  4,  a state  feedback  control  design  for 
a class  of  bilinear  systems  with  state-delay  is  presented.  The  stability  conditions 
derived  provide  a bound  on  the  system  delay,  and  define  an  attraction  region  based 
on  the  initial  condition.  The  future  work  proposed  for  consideration  is  presented  in 
Chapter  5. 


CHAPTER  2 

THE  NYQUIST  ROBUST  SENSITIVITY  MARGIN 
2.1  Introduction 

The  critical  direction  theory  introduced  by  Latchman  and  Crisalle  [35],  Latch- 
man  et  al.  [36],  and  later  generalized  by  Baab  et  al.  [2]  is  an  effective  approach 
for  analyzing  the  robust  stability  of  uncertain  systems  with  convex  and  non-convex 
uncertainty  value  sets.  A key  concept  introduced  by  the  theory  is  the  Nyquist  ro- 
bust stability  margin,  which  provides  a measure  of  robustness.  The  approach 

has  proven  useful  in  characterizing  the  robust  stability  of  single-input /single-output 
systems  with  real  affine  parametric  uncertainty,  among  others,  and  has  recently  been 
applied  to  the  design  of  robustly  stabilizing  i/oo  controllers  by  identifying  an  appro- 
priate weighting  function  for  the  controller  sensitivity  function  [29,  30]. 

This  chapter  proposes  an  alternative  robust-stability  analysis  which  has  the  ben- 
efit of  also  capturing  the  concept  of  robust  sensitivity,  hence  directly  incorporating 
the  notion  of  performance  robustness.  The  resulting  Nyquist  robust  sensitivity  mar- 
gin kN,s{oj)  is  inspired  on  the  critical  direction  theory  framework,  but  is  formulated 
to  take  into  account  in  an  explicit  fashion  the  effect  of  the  uncertain  systems  that 
have  the  worst-case  sensitivity. 

The  earlier  critical  direction  theory  involving  the  margin  k^iuj)  considers  only  a 
subset  of  the  uncertain  systems  in  the  robustness  analysis,  namely,  those  uncertain 
systems  whose  image  on  the  Nyquist  plane  lie  along  a pre-specified  oriented  line. 
Although  the  restricted  critical-set  of  systems  considered  leads  to  non-conservative 
conditions  for  robust  stability,  the  approach  ignores  all  perturbed  systems  that  have  a 
poor  sensitivity  (he.,  systems  located  close  to  the  critical  point  — 1-t-jO  on  the  Nyquist 
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plane)  whenever  these  lie  outside  the  oriented  line.  The  new  paradigm  involving  the 
margin  seeks  to  quantify  the  effect  of  the  systems  located  closest  to  the  critical 

point  through  the  introduction  of  a sensitivity  perturbation  radius  that  is  calculated 
at  each  frequency  by  solving  an  optimization  program. 

To  illustrate  the  approach,  the  robust  stability  analysis  proposed  is  developed 
for  uncertain  systems  described  by  rational  transfer  functions  with  real  affine  para- 
metric perturbations.  More  specifically,  the  numerator  and  denominator  polynomials 
depend  affinely  on  a set  of  real  parameters  that  are  known  to  belong  to  a given  uncer- 
tainty description.  A systematic  algorithm  for  the  calculation  of  kN,s{(^)  is  developed 
by  taking  advantage  of  simple  geometrical  features  adopted  by  the  Nyquist-plane 
images  of  such  systems  [21].  The  analysis  is  carried  out  in  detail  in  Section  4. 

The  robust  stability  of  the  real  affine  uncertain  systems  considered  in  Section  4 
can  be  analyzed  using  alternative  approaches,  for  example,  based  on  generalizations 
of  Kharitonov’s  methodology  [32].  In  particular,  one  may  adopt  the  approach  in 
Barmish  [3]  which  proposes  a strict  positivity  condition  that  must  be  evaluated  at 
a finite  number  of  frequencies,  or  the  box  theorem  [8],  or  the  worst-edge  algorithm 
of  Sideris  [48].  Furthermore,  for  the  robust  stability  of  interval  plants,  Wang  [53] 
has  shown  that  it  suffices  to  check  two  vertices.  Some  results  concerning  the  robust 
stability  of  control  systems  under  unstructured  as  well  as  parametric  uncertainty  have 
been  addressed  in  Chapellat  [9].  These  alternative  results  successfully  reveal  whether 
the  system  is  robustly  stable;  however,  in  contrast  to  the  Nyquist  robust  sensitivity 
margin  proposed  here,  they  do  not  provide  a scalar  indicator  of  the  closeness  to 
instability.  Hence,  the  scalar  k^s  can  be  used  to  compare  alternative  closed-loop 
designs  and  determine  a hierarchy  of  robust  stability  among  the  alternatives.  A 
recent  result  by  Wang  [54]  concerning  interval  plants  shows  that  the  maximum  i?oo 
norm  of  the  sensitivity  function  is  achieved  at  twelve  (out  of  sixteen)  Kharitonov 
vertices.  The  result,  however,  applies  to  interval  polynomials  while  our  approach 
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applies  to  transfer  functions.  A systematic  algorithm  for  the  calculation  of 
is  developed  by  taking  advantage  of  the  simple  geometrical  features  documented  in 
Fu  [21]  adopted  by  the  Nyquist-plane  images  of  such  systems. 

The  chapter  is  organized  as  follows.  In  Section  2,  the  classical  critical  direction 
theory  is  briefly  reviewed  for  contextual  reference.  Section  3 presents  the  definition  of 
the  new  Nyquist  robust  sensitivity  margin,  discusses  its  properties  and  computational 
challenges,  and  compares  and  contrasts  the  new  margin  with  its  Nyquist  robust 
stability  margin  predecessor.  The  application  of  the  Nyquist  robust  sensitivity  margin 
to  systems  with  affine  uncertainty  structure  is  presented  in  Section  4,  including  the 
details  of  a systematic  algorithm  for  the  efficient  calculation  of  the  margin.  Section  5 
presents  examples,  including  an  illustrative  case  showing  how  to  utilize  the  proposed 
method  for  calculating  a parametric  robust-stability  margin  that  is  interpreted  as  a 
blow-up  factor. 


2.2  Background 

A general  linear  time  invariant  system  (LTI)  can  be  represented  in  a state  space 
form  as  follows: 


x{t)  = Ax{t)  + Bu{t) 

y{t)  = Cx{t)  + Du{t)  (2.1) 

Furthermore,  the  representation  (2.1)  can  be  expressed  in  the  following  transfer  func- 
tion form: 

+ D = + D (2.2) 

provided  there  are  no  cancellations  between  the  numerator  and  denominator  poly- 
nomials. When  the  system  matrices  {A,  B,  C,  D)  are  uncertain,  the  transfer  function 
form  is  given,  for  the  MIMO  case,  by  G{s)  — Go('S)  + ^(^))  where  (?o('5)  is  a known 
transfer  matrix,  and  A(s)  is  the  transfer  matrix  representing  the  uncertainty.  The 
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transformation  of  system  (2.1)  into  (2.2)  allows  us  to  use  frequency  domain  tech- 
niques to  analyze  the  stability  and  performance  of  the  closed-loop  system.  Since  the 
development  of  the  Nyquist  robust  sensitivity  margin  in  this  chapter  requires  fre- 
quency domain  techniques  such  as  the  Nyquist  theorem,  the  transfer  function  form 
is  the  appropriate  environment  to  use  in  assessing  the  robust  stability  of  the  system. 

Consider  the  uncertain  single-input/single-output  transfer  function 

g{s)  = go{s)  + 6{s)  (2.3) 

shown  in  Figure  2.1,  where  go{s)  is  a nominal  system,  5(s)  6 A is  an  unknown 
perturbation  belonging  to  a known  set  of  allowable  perturbations  A.  The  closed- 
loop  system  of  Figure  2.1  is  said  to  be  robustly  stable  if  stability  is  ensured  for  all 
(i(s)  G A.  The  problem  under  consideration  is  the  analysis  of  the  robust  stabil- 
ity of  the  uncertain  closed-loop  system  (2.3)  under  negative  unity  feedback.  The 
developments  assume  the  following  standard  premises  that  are  commonly  used  in 
Nyquist-based  robustness  analysis:  (Al)  the  nominal  transfer  function  go{ju})  is  sta- 
ble under  negative- unity  feedback,  and  (^12)  the  uncertain  system  g{ju})  and  the 
nominal  system  go{jco)  have  the  same  number  of  open-loop  unstable  poles. 

-XO ► g(s)  ► 

f 


Figure  2.1:  The  uncertain  system  g{s)  = go{s)  -F  5{s)  in  a unity-feedback  configura- 
tion. 

The  key  concepts  and  definitions  pertaining  to  the  critical  direction  theory  are 
readily  summarized  utilizing  Figure  2.2.  First,  the  critical  line  is  the  oriented  line 
{i.e.,  a ray)  in  the  Nyquist  plane  originating  at  the  nominal  point  go{juj)  and  passing 
through  the  critical  point  —1  -f  jO.  The  critical  direction 
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(2.4) 


is  a unit-length  vector  with  origin  at  go{juj)  and  pointing  towards  the  critical  point. 
Then,  the  critical  ray  is  characterized  by  r{u)  = go{ju;)  +a  ddjuj)  for  a G . The 
uncertainty  value  set 


represents  the  Nyquist-plane  mapping  5f( jo;)  = go{j^)+d{ju})  of  the  uncertain  system. 
The  boundary  of  the  uncertainty  value  set  (2.5)  is  denoted  as  dV{u>).  Finally,  the 
critical  value  set  Vc{oj)  ■—  V(a^)  p]  r{co)  is  the  subset  of  V(u.))  that  lies  on  the  critical 
line. 

The  the  critical  value  set  Vc(u;)  may  be  convex,  i.e.,  a set  described  as  a single 
point  or  as  straight-line  segment  (such  as  the  straight-line  segment  go{juJi)gsU^i) 
joining  the  points  go{jiOi)  and  gsijoJi)  shown  is  Figure  2.2a),  or  nonconvex,  i.e.,  a 
union  of  isolated  points  and  straight-line  segments  (such  as  the  union  of  the  disjoint 
segments  go{j^^i)gi{j<-^i)  g2{j^i)g'i{j^i)  Figure  2.2b).  Note  that  it  is  possi- 

ble to  encounter  an  uncertain  system  with  a highly  nonconvex  value  set  V(w)  that 
nevertheless  features  a convex  critical  value  set  Vc(o;),  as  illustrated  in  Figure  2.2a. 

For  the  general  case  of  convex  or  nonconvex  critical  value  sets,  Baab  et  al.  [2] 


V(w)  = {g{juj)  : g{ju)  = goijuj)  + S{juj),  5(s)  G A} 


(2.5) 


define  the  critical  perturbation  radius 


(2.6) 


where 


(2.7) 
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d^UcOi)  Img(;<y) 


gs(M)  SiiM) 


Figure  2.2:  Uncertainty  value  sets  at  a frequency  uji.  (a)  convex  critical  value  set 
Vc{i^i),  (b)  nonconvex  critical  value  set  V{uJi).  Both  figures  show  the  worst-sensitivity 
plant  gsiji^i)-,  located  closest  to  the  point  —1  + jO. 

represents  the  minimal  distance  from  the  critical  point  —1  + jO  to  the  set  of  crit- 
ical boundary  intersections  Bc{u!)  :=  {5V(o;)  P|  ^0^^)}  \5o(^^),  where  '\'  is  the  set- 
exclusion  operator.  In  the  case  where  go{uj)  is  the  only  element  of  dV{u))f}r{juj), 
then  Bc{io)  :=  {po(^)}-  As  shown  in  Baab  et  al.  [2],  when  Vc{(^)  is  a convex  set,  as 
illustrated  in  Figure  2.2a,  the  definition  (2.6)  reduces  to 


pc{to)  ~ max{o;  G : g{juj)  = go{juj)  + a ddjco)  G Vc(o;)}  (2.8) 


a form  first  invoked  in  Latchman  et  al.  [36],  where  Pc(^)  is  simply  interpreted  as 
the  distance  between  the  critical  point  and  the  point  where  the  boundary  dV{u>) 
intercepts  the  critical  direction.  Finally,  the  Nyquist  robust  stability  margin  is  defined 
as 


kN{u) 


Pc{uj) 


\i  + go{ji^)\ 

The  main  result  of  Baab  et  al.  [2]  is  restated  in  the  following  theorem. 


(2.9) 


Theorem  1 Consider  the  uncertain  system  (2.3)  with  assumptions  (Al)  and  (A2). 
Then,  the  closed-loop  system  is  robustly  stable  under  unity  feedback  if  and  only  if 
ki^{uj)  < 1 Vo;. 
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Proof:  See  Baab  et  al.  [2].  B 

Note  that  the  theorem  is  valid  in  general  for  convex  as  well  as  nonconvex  critical 
value  sets  Vc(ci;).  Since  control  design  is  often  carried  out  under  sufficient-only  condi- 
tions, for  control  synthesis  purposes  it  may  be  acceptable  to  adopt  the  definition  (2.8) 
instead  of  (2.6)  when  working  with  nonconvex  critical  value  sets.  Then  the  resulting 
condition  k;^{uj)  < 1 Vw,  where  is  calculated  through  (2.9),  is  only  sufficient  for 
robust  stability. 


2.3  The  Nyquist  Robust  Sensitivity  Margin 
The  main  drawback  of  definition  (2.6)  is  that  the  resulting  Nyquist  robust  stabil- 
ity margin  value  kN{uj)  obtained  through  (2.6)  and  (2.9)  may  convey  no  information 
about  the  worst-case  sensitivity  in  the  value  set.  Figure  2.2b  shows  that  at  the  fre- 
quency u)  = uji  the  plant  gs{j^i)  is  the  element  of  V{uJi)  that  is  closest  to  the  point 
-I  + jO.  Hence  the  sensitivity  is  the  largest  among  all  the  plants  in  the 

value  set.  Note  that  since  gsijuJi)  ^ Vc(wj),  this  plant  is  ignored  in  the  classical 
critical  direction  analysis  presented  in  Section  2 which  focuses  only  on  plants  that  lie 
along  the  critical  direction. 

In  this  section  an  alternative  approach  is  presented  to  include  sensitivity  effects 
in  the  robustness  margin.  To  this  end  we  define  the  sensitivity  perturbation  radius 


|1 + 5o(j^^)l  - ^(t^)  if  -l+jO^V{u) 
11  + go{j‘^)\  + otherwise 


(2.10) 


where 


niw)  = min  |l  + -z|  (2-11) 

zedV{<.) 

represents  the  minimum  distance  between  the  critical  point  — l-fjO  and  the  boundary 
set  dV{w).  Then,  in  a fashion  analogous  to  (2.9),  the  Nyquist  robust  sensitivity  margin 
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is  defined  as 


(2.12) 


Img(70) 


Reg(jco) 


Figure  2.3:  Illustration  of  the  inverse-sensitivity  circle  of  radius  t}{uj)  introduced  in 
definition  (2.11). 

Figure  2.3  gives  an  interpretation  of  r]{u)  defined  in  (2.11)  as  the  radius  of 
the  inverse-sensitivity  circle,  namely,  the  smallest  circle  with  center  at  -1  -\-  jO  that 
contains  a point  belonging  to  the  boundary  dV{io).  Furthermore,  the  definition  (2.11) 
and  Figure  2.3  can  be  used  to  conclude  that  r]{oj)  = |1  + 5s(j^)|)  where  Qsij^)  is  the 
perturbation  in  V(o;)  that  has  the  worst  sensitivity. 

It  is  also  of  interest  to  note  that  = 1 corresponds  to  the  case  where 

—1  -\-j0  G dV{Lo).  This  follows  from  the  fact  that  fcjv,s(n;)  = 1 if  and  only  if  r]{oj)  = 0, 
and  the  latter  equality  is  realized  from  the  optimization  problem  (2.11)  only  when 
-1  -h  jO  G dV{uj).  Finally,  it  is  of  utility  for  the  suite  to  note  that  at  all  frequencies  oo 


This  inequality  is  derived  as  follows.  Since  the  value-set  boundary  dV{u))  contains  as 
a subset  the  set  of  critical  boundary  intersections  Bc{io),  the  optimization  problem 
(2.11)  is  carried  out  over  an  optimization  domain  that  is  a superset  of  the  optimization 


T]{uj)  < ^{uj) 


(2.13) 
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domain  used  in  the  optimization  problem  (2.7).  Consequently,  the  solutions  to  the 
respective  optimization  problems  must  follow  the  relationship  (2.13). 

Theorem  2 Consider  the  uncertain  system  (2.3)  with  assumptions  (Al)  and  (A2). 
Then,  the  closed-loop  system  is  robustly  stable  under  unity  feedback  if  and  only  if 
kN,s{^)  < 1 

Proof:  From  the  zero-exclusion  principle  [4]  it  can  be  claimed  that  the  uncertain 
system  (2.3)  under  assumptions  (Al)  and  (A2)  is  robustly  stable  if  and  only  if  —1  + 
jO  ^ V(a»).  Therefore,  it  must  be  shown  that  under  the  definitions  (2.10)-(2.12) 
for  ps{to),  the  condition  < 1 is  equivalent  to  the  set  membership  condition 

-1+jO^  V(a;). 

First,  to  prove  sufficiency  one  must  show  that  kN,s{'^)  < 1 Vw  implies  that 
—1  -I-  jO  ^ V(o;).  The  proof  proceeds  by  contradiction.  Assume  that  fc;v,s(‘^)  < 1 and 
that  there  exists  a frequency  w such  that  —1  + jO  G V(o;).  Invoking  the  sensitivity- 
perturbation  radius  expression  (2.10)  for  the  case  where  —1  + jO  G V{lo)  and  the 
definition  (2.12)  it  follows  that 


kN,s{pj) 


_ 


|i+go(jt^)l+y?(t*^) 

|l+9o(jo))| 


1 + 


ll+9o(jw)|  11+90(11^)1  ' |l+9o(li^)| 

Since  by  definition  77(0;)  > 0,  then  the  equation  above  implies  that  kM,s{<^)  > 1,  which 
is  a contradiction.  This  proves  sufficiency. 

Second,  to  prove  necessity  one  must  show  that  at  any  frequency  u)  the  condition 
-1+jO  i V(w)  implies  that  fc^,,(a;)  < 1.  Assume  that  -1+jO  ^ V(o;).  Invoking  the 
sensitivity-perturbation  radius  expression  (2.10),  now  for  the  case  where  — 1 -h  jO  ^ 
V(o;),  and  the  definition  (2.12),  it  follows  that 


kN,s{<^)  — lY 


Ps{i^) 


+ 9oU^)\ 


|1  + yo(ja;)|  -77(07) 

|1  + 5o(iw)l 


77(0;) 


11  + 5o(i^^)| 


(2.14) 


Since  in  this  case  -1  + jO  ^ V(o;),  it  follows  that  -1  + jO  ^ dV{u),  and  hence  from 
(2.11)  it  is  concluded  that  t]{uj)  > 0.  Furthermore,  from  (2.11)  it  is  obvious  that 
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r]{uj)  < |1  + 5o(w)|.  Hence,  it  follows  that  0 < H which  can  be  used  in 

to  conclude  that /cjv,s(<^)  < 1-  * 

Figure  2.2a  illustrates  a special  situation  where  /cAr,s(o;)  = /cat (a;).  This  follows 
from  a simple  argument  using  the  elements  shown  in  the  figure,  where  is  clear  that 
in  this  case  t]{u)  — ^{to)  = |1  + gs{j^)\-  Hence  Ps{co)  = Pc(^)  from  (2.10)  and 
(2.6),  and  therefore  it  follows  from  (2.12)  and  (2.9)  that  fcAr,s(ct))  = k!^{uj).  It  is 
also  straightforward  to  verify  that  the  robustness  margins  satisfy  the  following  two 
properties:  (PI)  if  kN^s{<^)  < 1)  then  /cjv,s(^)  ^ ^jv(w),  and  (P2)  if  kN^s{<^)  > then 
kN,s{^)  — ^iv(t^)-  These  two  properties  follow  in  a straightforward  fashion  after  using 
the  inequality  (2.13),  the  perturbation  radius  definitions  (2.10)  and  (2.6),  and  the 
robustness  margin  definitions  (2.12)  and  (2.9)  . Although  in  general  kN,s{‘^)  7^  kN{uj), 
as  suggested  in  Figure  2.2b,  both  margins  are  nevertheless  equivalent  as  indicated  in 
the  following  theorem. 

Theorem  3 The  Nyquist  robust  sensitivity  margin  kN,s{‘^)  the  Nyquist  robust 
stability  margin  A)Ar(w)  are  equivalent  in  the  sense  that  (i)  fciv,s(<^)  < 1 kN{uj)  < 
1,  (a)  kN,s{i^)  — 1 kN{(^)  = I,  and  (Hi)  kN,s{oj)  > 1 k^ita)  > 1- 

Proof:  The  proof  of  sufficiency  is  developed  below  for  cases  (i)-(iii).  The  proof 
of  necessity  for  the  three  cases  in  question  follows  an  analogous  argument,  and  is 
therefore  omitted  here  for  brevity.  For  case  (i),  assume  kf^^s{oj)  < 1 and  utilize  (2.12) 
to  conclude  that 

Ps{uj)  < \1  + go{jaj)\  (2.15) 

Also,  from  Theorem  2 and  from  the  zero-exclusion  principle  [4],  the  condition  k^^sioj)  < 
1 implies  that  — 1-l-jO  ^ V(a;);  hence  from  equations  (2.6)  and  (2.10)  the  appropriate 
expressions  for  the  respective  perturbation  radii  are  Pc{oj)  = |1  + go{joj)  \ — ^(w)  and 
Ps{tv)  = |1  -I-  goU^)  \ ~ From  the  latter  two  equations  and  inequality  (2.13)  it 

follows  that 


Pc{0j)  < Ps{0}) 


(2.16) 
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Inequalities  (2.15)  and  (2.16)  imply  that  pdco)  < |1  + 9o{j‘^)\  which  yields  the  result 
kM{u})  < 1 after  invoking  (2.9).  For  case  (ii),  assume  — 1 and  utilize  (2.12) 

to  conclude  that  Ps(w)  = |1  + Po(^)|,  which  in  turn  from  equation  (2.10)  implies  that 
T]{u;)  = 0.  Since  ry(a;)  = 0 solves  the  optimization  problem  (2.11),  it  follows  that 
-1+jO  G dV{uj).  Now,  the  fact  that  -1+jO  G dV{ij)  implies  that  -1+jO  G Bc{oo), 
where  Bc{co)  is  the  optimization  domain  in  (2.7).  Since  —1  + jO  G Bc{uj),  it  follows 
that  the  solution  to  the  optimization  problem  (2.7)  is  ^(w)  = 0,  which  can  be  used 
to  conclude  from  (2.6)  that  pc{uj)  = |1  + Po(<^)|-  Substitution  the  latter  equality 
into  (2.9)  yields  kN{uj)  = 1.  For  case  (iii),  assume  A:;v,s(<^)  > 1 and  utilize  (2.12)  to 
conclude  that 

Ps{uj)  > \1  + go{ju})\  (2.17) 

Also,  from  Theorem  2 and  from  the  zero-exclusion  principle  [4],  the  condition  kpf^si^)  > 
1 implies  that  -l+jO  G V(o;);  hence  from  equations  (2.6)  and  (2.10)  the  appropriate 
expressions  for  the  respective  perturbation  radii  are  Pc{oj)  = |1  + 9o{j^)  \ +^(^)  and 
Ps{uj)  = |1  + 9o{ji^)  \ + v{^)-  From  the  latter  two  equations  and  inequality  (2.13)  it 
follows  that 

Pc{i^)  > Ps{i^)  (2-18) 

Inequalities  (2.17)  and  (2.18)  imply  that  Pc{uj)  > |1  -I- po(i<^)|  which  yields  the  result 
k^ii^)  > 1 after  invoking  (2.9).  ■ 

The  Nyquist  robust  sensitivity  margin  serves  a role  analogous  to  that  of 

the  structured  singular  value  p{co)  [17]  or  the  multivariable  stability  margin  km{i^) 
[46],  as  a scalar  indicator  of  robust  stability.  Given  that  the  optimization  problem 
(2.11)  must  be  solved,  the  deployment  of  an  analysis  approach  based  on  /cw,s(^^) 
requires  knowledge  the  value-set  boundary  dV{uj).  Fortunately  this  information  is 
available  in  a number  of  problems  of  interest,  such  as  the  case  of  systems  with  real 
affine  uncertainty  structure  discussed  in  the  following  section. 
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2.4  Application  to  Systems  with  Affine  Uncertainty  Structure 
The  robust  stability  analysis  approach  proposed  is  applied  to  a class  of  uncertain 
systems  with  real  affine  uncertainty  structure  of  the  form 


_ no{s)  + Y7i=iQMs) 

do{s)  + EUQ^d^{s) 


(2.19) 


where  rii{s)  are  numerator  polynomials  of  known  order  i and  known  real  coefficients 
Uik,  k = 0, i = 0,1,..., p,  and  where  d,(s)  are  denominator  polynomials  of 
known  order  m and  known  real  coefficients  dik,  k = 0,1,..., m,  i — 0,  l,...,p.  The 
element  q G Q is  a vector  of  real  perturbation  parameters,  where  the  real  uncertainty 
domain 

Q = {q  € 7^P  : qf  < q < q,^,  i = 1,  2, . . .,  p}  (2.20) 


is  a bounded  rectangular  polytope.  In  this  case  the  uncertainty  value  set  V(o;)  is 
simply  the  map  g{juj,  Q)  : TV'  x Q ^ C. 

The  objective  is  to  calculate  the  value  of  /cAr,s(w)  as  a function  of  frequency 
using  the  expression  (2.12).  This  in  turn  requires  the  calculation  of  the  sensitivity 
perturbation  radius  Ps{oo)  through  its  defining  equation  (2.10).  Note  that  in  order 
to  apply  (2.10)  two  problems  must  be  addressed,  namely,  the  optimization  program 
(2.11)  must  be  solved  to  find  the  inverse-sensitivity  radius  r]{u)  {Problem  I),  and  the 
set-membership  clause  -1  + jO  ^ V(u;)  must  be  assessed  as  true  or  false  {Problem 
II)  so  that  the  appropriate  branch  of  equation  (2.10)  can  be  identified. 

It  is  shown  in  [21]  that  the  mapping  g{joj,  E{Q)  (denoted  in  the  suite  as  the 
value-set  frame  at  the  frequency  lo),  where  E{Q)  represents  the  set  of  edges  of  Q, 
spans  the  boundary  set  dV{oo).  Furthermore,  the  frame  g{juj,  E{Q))  is  a set  comprised 
of  arcs  of  circles  and  straight-line  segments  [21].  More  precisely,  let  Ei{Q)  and  its 
corresponding  extreme  points  q~  and  qf  represent  the  i — th  edge  of  the  rectangular 
polytope  Q.  Then  the  frame  g{ju),E{Q))  is  composed  of  a set  of  frame-elements 
g{ju),Ei{Q)),  and  each  frame  element  is  either  a straight-line  segment  or  an  arc 
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of  a circle.  These  simple  geometric  properties  of  the  frame  allow  the  development 
of  a precise  solution  of  Problem  I.  In  fact,  the  minimization  problem  (2.11),  which 
is  equivalent  to  finding  the  minimum  distance  between  the  point  —1  + jO  and  the 
boundary  of  the  value  set,  reduces  to  a simple  geometric  problem:  finding  the  shortest 
distance  between  the  point  -1  + jO  and  an  arc  of  circle  or  a straight-line  segment. 
Problem  (2.11)  can  then  be  posed  for  each  frame  element,  and  the  smallest  solution 
found  after  considering  all  the  edges  of  Q yields  the  value  of  77(0;)  sought. 

For  completeness  it  is  convenient  to  briefly  summarize  relevant  geometrical  con- 
cepts regarding  lines  and  arcs  of  circles.  The  interested  reader  is  referred  to  [15]  for 
further  details.  A line  passing  through  points  pi,p2  ^ C is  defined  by 

L{Pi,P2)  :={z  eC  : z = pi  +u{p2  -pi),w  G 77} 


and  a circle  with  radius  r and  center  zq  is  given  by 


C{r,  Zq)  :=  {2:  G C : |z  - 2:qP  == 


The  arc  A{p\,p2,pz)  of  a supporting  circle  C{r,Zo)  is  described  by  three  points. 

One  important  issue  to  resolve  is  whether  the  map  g{juj,  Ei{Q))  of  a given  edge 
Ei{Q)  is  a straight-line  or  an  arc  of  a circle.  This  can  be  resolved  by  taking  advantage 
of  the  cross  product 


Pi  X P2 


Re{pi)  Rc{p2) 

Im{pi)  Im{p2) 

where  | • | represents  the  determinant  operator.  Selecting  three  distinct  points  Pi,P2, 
and  p3  of  the  map  g{jco,  Ei{Q)),  if  follows  that  if  {p^-pi)  x {p2~Pi)  < 0 (>  0)  it  can  be 
concluded  that  the  segment  is  an  arc  of  a circle  turning  to  the  right  (left).  If  the  cross 
product  equals  zero,  then  the  three  points  are  collinear  and  the  segment  is  a straight 
line.  Two  of  the  points  in  question  should  be  pi  — g{jco,  q~)  and  p^  = g{jco,  q/)-  The 
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?2 


Figure  2.4:  The  center  zq  and  radius  r of  the  supporting  circle  of  the  arc  A(pi,p2,P3) 
are  determined  from  the  intersection  of  the  auxiliary  lines  Li  and  L2 

third  point  can  be  taken  as  the  image  of  a distinct  point  on  the  edge  Bi(Q)  that  can 
be  selected  arbitrarily,  say  for  example,  p2  ~ g{juj,  {q^  + q^)/2). 

For  the  case  where  the  map  g{ju),  Ei{Q))  of  an  edge  Ei{Q)  is  an  arc,  the  minimum 
distance  from  —1  + jO  to  the  arc  can  be  given  by  either  (z)  the  distance  to  one  of  the 
end  points  of  the  arc,  or  (ii)  the  distance  to  an  internal  point  of  the  arc.  Clearly,  if  the 
ray  originating  at  the  center  of  the  supporting  circle  of  the  arc  and  passing  through 
— 1 + jO  does  not  intersect  the  arc,  then  the  minimum  distance  can  be  determined 
from  one  of  the  two  end  points  of  the  arc.  On  the  other  hand,  if  the  ray  intersects 
the  arc,  then  the  distance  between  — 1 + jO  and  the  point  of  intersection  defines  the 
minimum  distance  sought. 

Finally,  the  procedure  described  requires  finding  the  supporting  circle  of  an  arc. 
From  Figure  2.4,  the  center  zq  of  the  supporting  circle  of  an  arc  that  passes  through 
three  distinct  points  pi,P2,  and  ps  on  the  complex  plane  can  be  determined  from  the 
intersection  of  the  lines  Li  :=  (pi  + P2)/2  + jui{p2  — Pi),  and  L2  :=  (ps  + P2)/2  + 
ju2{pz  — P2),  where  U2  G TZ.  The  radius  r is  then  found  in  an  obvious  fashion,  say 
for  example  as  r = \p2  — zq\. 

For  the  case  where  the  map  g{jui,  Ei{Q))  of  an  edge  Ei{Q)  is  a straight-line 
segment,  the  minimum  distance  to  the  point  —1  + jO  is  found  using  a procedure 
formally  similar  to  the  case  of  the  arcs.  First  a supporting  line  is  found.  Then,  one 
finds  the  point  of  intersection  between  the  supporting  line  and  a normal  line  that 
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passes  through  — 1+jO.  The  intersection  point  gives  the  minimal  distance  to  — 1+jO 
if  the  intersection  point  is  also  an  element  of  the  straight-line  segment.  Otherwise, 
the  minimum  distance  is  defined  as  the  distance  between  — 1 -I- jO  and  one  of  the  two 
end-points  of  the  straight-line  segment. 

The  procedure  described  above  solves  Problem  I,  yielding  a numerical  value  for 
the  inverse-sensitivity  radius  t]{uj)  at  each  frequency.  Problem  II  can  be  addressed 
efficiently  through  the  assistance  of  the  following  theorem,  which  is  a restatement  of 
an  equivalent  theorem  derived  in  Baab  et  al.  [2].  A detailed  proof  is  given  in  the 
original  reference. 

Theorem  4 Consider  the  real-affine  uncertain  system  (2.19)-  (2.20)  configured  in 
the  unity-feedback  form  given  in  Figure  2.3  under  the  assumptions  (Al)  and  (A2). 
Then  -1  jO  ^ V(o;)  if  and  only  if  at  frequency  u the  following  linear  equal- 
ity/inequality problem  is  infeasible: 


Aq  = b 


(2.21) 


subject  to 


Bq<  b+ 


(2.22) 


where 
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where  [■]  represents  the  greatest-integer  function. 

Proof:  See  Baab  et  al.  [2],  B 

In  summary,  for  the  system  (2.19)  with  parametric  uncertainty  (2.20)  it  is  possi- 
ble to  solve  Problem  I and  calculate  with  very  high  numerical  precision  the  sensitivity 
radius  rj{uj)  because  the  solution  to  (2.11)  is  given  by  a set  of  simple  algebraic  equa- 
tions. In  addition,  it  is  possible  to  solve  Problem  II  in  a numerically  efficient  fashion 
because  the  condition  -1  -I-  jO  G V(w)  can  be  determined  via  a simple  feasibility 
problem  involving  linear  equalities  and  inequalities.  Hence,  the  sensitivity  pertur- 
bation radius  ps{^)  in  (2.10)  and  the  Nyquist  robust  sensitivity  margin  fc;v,s(a;)  in 
(2.12)  can  be  computed  precisely  and  efficiently. 


Three  examples  are  presented.  The  first  one  calculates  the  margins 
and  and  to  compare  and  contrast  their  values  and  to  shed  light  into  their 

interpretation.  The  second  example  illustrates  in  a dramatic  fashion  the  fact  that 
kN,s{^)  provides  a more  meaningful  indication  of  the  degree  of  robust  sensitivity  of 
the  closed  loop.  Finally,  the  last  example  is  designed  to  illustrate  how  the  concepts 
proposed  here  can  be  utilized  to  formulate  and  calculate  an  alternative  robustness 
measure,  namely  a parametric  stability  margin. 

2.5.1  Example  1 

Consider  the  affine  system  of  the  form  (2.19)  with  the  structure  [29] 


2.5  Examples 


g{s,q)  c(s) 


5 s -h  Qi 


(2.23) 


s^-hQ2  s-hqs 
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where 

3603.7935  s + 18018.9673 
“ fi2  + 1434.5016  s - 2312.4499 

is  a feedback  controller.  Let  the  real  perturbation  parameters  belong  to  the  uncer- 
tainty domain 

Q = {(71,72,93)  G 77^  : 0 < < 8,  -2  < 92  < 6,  -19  < 93  < -11}  (2.24) 

Figure  2.5  shows  the  uncertainty  value  set  for  (2.23)  at  the  frequency  a;  = 9,  including 
the  corresponding  sensitivity  circle  centered  at  —1  -I-  jO.  Note  that  Figure  2.5  also 
shows  the  frame  of  the  value  set  of  (2.23),  namely,  the  straight-line  segments  and 
arcs  of  circles  that  result  from  the  mapping  of  all  the  edges  of  Q.  The  problem  is 
to  analyze  the  robust  stability  of  the  feedback  loop  involving  the  uncertain  system 
(2.23)  subject  to  the  uncertainty  description  (2.24). 

The  margin  fciv,s(w),  calculated  following  the  algorithm  given  in  Section  4,  and 
the  margin  /cat (a;),  calculated  using  the  technique  described  in  Baab  et  al.  [2],  are 
plotted  in  Figure  2.6  for  frequencies  lo  G [10~^,  10].  Given  that  fcAr,s(a;)  < 1 Vo;, 
it  can  be  concluded  from  Theorem  2 that  the  closed-loop  system  is  robustly  stable. 
Since  the  two  margins  are  equivalent,  the  values  of  < 1 reported  in  the  figure 

correspond  to  values  kN{u)  < 1 at  the  same  frequency,  consistent  with  Theorem  3. 
Note  that  Figure  2.6  shows  that  in  this  particular  case  kN,s{<^)  is 

an  upper  observation  is  consistent  with  property  (PI)  which  implies  that  is 

an  upper  bound  for  A:jv(a>)  when  the  system  is  robustly  stable. 

2.5.2  Example  2 

Consider  the  system  [2] 

9(s,g)=c(s)^  (2.25) 

where  n{s)  = -f  (4  -I-  O.491  -1-  0.292)5  -I-  (20  -f  91  —93),  d{s)  = -I-  (9.5  -I-  O.591  — 

0.592  + 0.593)5^  -h  (27  +291  -b  92)5^  -b  (22.5  - 9i  -b  93)5  + 0.1,  c{s)  = 0.3s  -b  1,  and  where 
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Figure  2.5:  Frame  for  the  value  set  of  system  (2.23)  at  cj  = 9,  and  the  corresponding 
inverse-sensitivity  circle.  The  nominal  plant  is  indicated  by  the  ’-f  ’ marker. 


Figure  2.6:  Value  of  and  as  a function  of  frequency  for  the  first  example. 
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the  real  perturbation  parameters  belong  to  the  polytope 

Qi  = {(91,92,93)  en^:-3<qi<3,i  = 1,2,3}  (2.26) 

Introducing  a parametric  blow-up  factor  a = 1.86  to  define  a larger  uncertainty 
domain 

Q = a Qi  (2-27) 

such  that 

Q = {(91, 92, 93)  e : -5.58  <qi<  5.58,  z = 1,  2,  3} 

and  using  lo  - 4.72  yields  the  frame  g{juj,E{Q))  depicted  in  Figure  2.7.  The  re- 
sults obtained  are  as  follows:  Ps{uj)  = 0.8466,  Pc(^^)  = 0.0206,  A: tv,s(<^)  = 0.8577,  and 
kN{uj)  = 0.0209.  Both  /cAr,s(a;)  and  are  less  than  unity  at  the  frequency  con- 

sidered. Further  analysis  shows  that  the  result  holds  at  all  frequencies;  hence,  the 
closed  loop  system  is  robustly  stable  according  to  Theorem  2.  Certainly,  kN,s{^^) 
and  kN{uj)  are  equivalent  in  terms  of  assessing  the  robust  stability  of  the  system,  as 
proved  in  Theorem  3.  However  in  the  terms  of  robust  performance,  A:Ar,s(a;)  is  a more 
meaningful  metric  than  kN{uj)  because  the  larger  value  of  k^^si^)  better  reflects  the 
fact  that  the  value  set  in  Figure  2.7  is  in  close  proximity  of  the  critical  point  -1  + jO. 
In  fact,  at  the  frequency  in  question  k^^si,^)  — 0.8577  is  a much  larger  value  than 
k!^{uj)  = 0.0209.  The  significantly  smaller  value  of  k^iuj)  does  not  yield  compa- 
rable insight  into  the  proximity  of  the  worst-sensitivity  perturbation  to  the  point 
-1  -h  jO  because  the  calculation  of  does  deliberately  ignores  all  uncertainties 

lying  outside  the  critical  direction. 

2.5.3  Example  3 

An  useful  application  of  the  Nyquist  robust  sensitivity  margin  proposed  is  the 
calculation  of  a parametric  stability  margin.  Consider  the  system  (2.25)  used  in 
Example  2 and  the  polytope  (2.26)-(2.27)  featuring  a variable  blow-up  factor  a;  > 0. 


25 


Figure  2.7:  Frame  for  the  value  set  of  system  (2.25)  at  a;  = 4.72  and  a — 1.86.  The 
nominal  plant  go{juj)  is  indicated  by  the  ’+’  marker. 

Figure  2.8  shows  the  Nyquist  robust  sensitivity  margin  kM,s  '■=  maxkM,s{^)  that 

UJ 

results  when  magnifying  the  original  perturbation  polytope  Q by  different  blow-up 
factor  values.  The  numerical  study  shows  that  when  the  blow-up  factor  has  the  value 
a = 1.89  then  the  Nyquist  robust  sensitivity  margin  kN,s  is  approximately  equal  to 
unity,  hence  reaching  the  limit  of  robust  stability.  The  limiting  value  a = 1.89  is  the 
parametric  robust  stability  margin  for  the  uncertain  closed  loop.  In  other  words,  the 
controller  c(s)  introduced  in  Example  2 can  robustly  stabilize  the  closed-loop  system 
subject  to  a parametric  blow-up  factor  of  the  parametric  uncertainty  domain  (2.27) 
less  than  a.  Note  that  the  blow-up  factor  used  in  Example  2 is  o;  = 1.86  < a;  hence, 
the  uncertain  closed-loop  of  Example  2 is  robustly  stable. 
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Figure  2.8:  Plot  of  the  Nyquist  robust  sensitivity  margin  = maxkM,s{(^)  as  a 
function  of  the  blow-up  factor  a.  The  parametric  robust  stability  margin  is  a = 1.89, 
which  corresponds  to  the  value  of  the  blow-up  factor  a that  makes  k^^g  approximately 
equal  to  unity. 

2.6  Conclusions 

The  new  concept  of  a Nyquist  robust  sensitivity  margin  can  be  used  to  quantify 
the  robust  stability  of  uncertain  closed-loop  systems  while  at  the  same  time  produc- 
ing a meaningful  indication  of  the  worst-case  sensitivity  that  is  realized.  Hence,  in 
this  sense  the  approach  is  more  attractive  than  the  classical  Nyquist  robust  stability 
margin  framework,  which  ignores  all  systems  that  do  not  lie  along  the  critical  direc- 
tion, and  that  may  therefore  exclude  from  the  analysis  perturbed  systems  that  have 
the  worst  sensitivity. 

On  the  other  hand,  in  general  the  calculation  of  the  Nyquist  robust  sensitivity 
margin  may  involve  more  numerically  intensive  optimization  work,  since  the  program 
(2.11)  is  a superset  of  the  program  (2.7).  In  other  words,  the  calculation  of  kN^g{uj) 
requires  knowledge  of  the  entire  value-set  boundary  5V(u;),  whereas  the  calculation 
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of  kN{u})  requires  knowledge  of  only  those  points  of  dV{u)  that  lie  along  the  critical 
direction. 

The  examples  presented  illustrate  the  ability  of  the  Nyquist  robust  sensitivity 
margin  methodology  to  produce  meaningful  quantitative  measures  of  robustness  for 
uncertain  systems,  even  in  the  case  where  the  uncertainty  value  set  is  originated 
by  a real  parametric  uncertainty  description.  The  examples  also  shows  how  the 
proposed  paradigm  can  be  used  to  characterize  alternative  robustness  measures,  such 
as  a parametric  blow-up  factor  for  a real  uncertainty  description  comprised  of  a 
rectangular  polytope. 

The  numerical  algorithm  used  to  solve  the  problem  of  Section  4 calls  for  a mod- 
est computational  requirement  because  the  calculation  of  the  sensitivity  radius  can 
be  carried  out  in  a straightforward  fashion.  It  is  anticipated,  however,  that  the 
computational  cost  associated  with  other  particular  real  parametric  uncertainty  de- 
scriptions may  become  significantly  more  expensive  given  that  the  general  parametric 
uncertainty  analysis  problem  is  found  to  be  NP  hard  [44]. 

2.7  Supplementary  Calculation  Algorithms 

Further  detail  on  the  computational  techniques  discussed  in  Section  2.4  are  pre- 
sented in  this  section.  First,  an  alternative  algorithm  to  finding  the  supporting  circle 
of  an  arc  is  introduced.  Second,  a simple  algorithm  to  find  the  minimum  distance 
between  the  critical  point  and  a line  is  discussed.  Finally,  a technique  is  discussed 
for  determining  whether  the  intersection  point  of  the  line  through  the  critical  point 
and  the  supporting  circle  of  the  arc  segment  actually  lies  on  the  arc. 

2.7.1  Supporting  Circle  of  an  Arc 

An  alternative  method  for  finding  the  supporting  circle  on  an  Arc  defined  by 
three  points,  A(rri,  X2,  X3),  is  to  utilize  the  equation  of  the  circle.  A circle  that  passes 
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through  the  point  x and  centered  at  the  point  C is  given  by 

{Xr  - Crf  + {Xi  - CiY  = (2.28) 

where  the  subscript  r refers  to  the  real  part  of  the  complex  point  x,  and  the  subscript 
i refers  to  its  imaginary  part.  Now,  using  the  three  points  that  define  the  Arc,  namely 
Xi,X2,  and  X3,  into  equation  (2.28),  the  following  three  equations  are  obtained: 

— ‘lXy\Cf  + — 2xiiCj  -)- 

X^2  ~ 2Xr2Cr  + X^2  — 2Xi2Cj  + Cf  = 

X^3  — 2XrzCr  + + X^3  — 2Xi^C{  + Cf  = (2.29) 

The  three  equations  (2.29)  have  three  unknowns,  namely  C'r,Cj,  and  r,  which  pre- 
cisely identify  the  supporting  circle  of  the  Arc  defined  by  Xi,X2,  and  X3. 

2.7.2  Minimum  Distance  between  a Line  and  a Point 

Given  that  uncertain  systems  of  the  form  (2.19)  produce  either  a line-segment 
or  an  Arc,  it  is  important  to  be  able  to  find  the  minimum  distance  between  the 
critical  point  —1  + jO  and  a line-segment.  The  projection  technique  can  be  utilized 
to  do  just  that.  The  following  steps  describe  the  procedure:  First,  given  the  critical 

point  cp,  and  the  two  end-points  ps  and  pe  of  the  line-segment,  three  vectors  are 

defined  as  follows: 


Re{cp) 

Re{ps) 

Re{pe) 

^cp  — 

Im{cp) 

Vps  — 

Im{ps) 

'^pe  — 

Im{pe) 

Next,  define  the  vector  v = Vpe  — Vps,  and  the  the  direction  vector  Vd  — pe  — ps. 
Finally,  the  projection  is  computed  as 

p {pcp  '^ps)  ■ ^ 
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where  the  refers  to  the  dot  product  operator.  Now,  the  minimum  distance  is 
calculated  based  on  the  sign  of  the  projection  as  follows: 

- if  P < 0 then  dmin  = \cp  — ps\ 

- if  P > 1 then  dmin  = \cp  — pe\ 

- if  0 > P < 1 then  d^in  = \cp  - {ps  + proj{pe  - ps))| 

2.7.3  Identifying  Points  on  the  Arc 

One  decision  that  has  to  be  made  in  Section  2.4  is  whether  the  intersection  point 
of  the  ray  originating  from  the  center  of  the  supporting  circle  of  the  Arc  and  passing 
through  the  critical  point  -1  + jO  lies  on  the  Arc.  The  key  step  is  to  identify  the 
valid  Arc  phase  range,  /arc,  such  that  if  a point  p G larc  then  p lies  on  the  Arc.  The 
algorithm  can  be  summarized  as  follows: 

• Shift  the  supporting  circle  and  the  intersection  point  p to  the  origin. 

• Find  the  phase  of  the  start  and  end  points  of  the  Arc,  and  convert  them  such 
that  the  phase  ranges  between  [0  , 27t]. 

• Denote  the  smallest  phase  (j)min  and  the  largest  <t>max- 

• Find  the  phase  of  a mid  point  on  the  Arc  and  denote  it  (j)m- 

- If  (prnin  ^ 0m  ^ 0mai  then  larc  ~ [0min  ) 0max]-  Furthermore,  if  j?  G larc 

then  p lies  on  the  Arc. 

- Else  if  (f>nfi  ^ 0mm)  and  (j)jji  ^ 0mai)  then  larc  [0  ) 0mm]  U[0mai  ) Oj- 
Thus,  if  p G larc  then  p lies  on  the  Arc. 

- Else  larc  = [0mai  , 4>m]  U[0m  , 0]  |J[0  , 0mm]-  Therefore,  if  p G larc  then 
p lies  on  the  Arc. 


CHAPTER  3 

SLIDING  MODE  CONTROL  FOR  TIME-DELAY  SYSTEMS 


3.1  Introduction 

Delay  is  inherit  in  some  control  systems  such  as  processes  involving  heat  or  mass 
transport.  The  presence  of  delay  in  a dynamic  system  can  have  a destabilizing  effect 
or  can  cause  poor  performance.  Furthermore,  delay  can  pose  a significant  challenge 
to  ensure  closed  loop  stability  [27].  Throughout  the  literature,  a variety  of  linear  and 
nonlinear  controllers  have  been  used  to  stabilize  time-delay  systems,  where  the  delay 
may  appear  in  the  state,  input,  or  in  both.  Local  and  global  stability  conditions  have 
also  been  derived  to  guarantee  the  asymptotic  stability  of  the  closed  loop  system. 

The  emphases  of  this  work  in  on  sliding  mode  control  (SMC),  a technique  known 
for  its  robustness  with  respect  to  perturbations  and  system  uncertainties,  that  has 
been  used  to  stabilize  systems  with  time  delays;  however,  most  of  the  literature 
focuses  on  systems  with  either  state  delay  [6,  33,  43,  28]  or  with  input  delay  [25,  45]. 
Some  work  has  been  done  regarding  systems  with  simultaneous  state  and  input  delays 
[22,  20,  55].  In  Gouaisbaut  et  al.  [22]  a sliding  mode  controller  is  designed  to  stabilize 
a linear  systems  with  input  and  state  delay.  The  technique  is  based  on  transforming 
the  system  into  a regular  form  [31],  and  then  a memory  control  law  that  depends  on 
previous  values  of  the  input  is  designed  the  ensure  reaching  the  manifold  as  well  as 
the  asymptotic  stability  of  the  closed-loop  system.  Lyapunov-Kharsovski  methods 
are  employed  to  derived  the  stability  conditions.  The  method  presented  in  Feiqi  et 
al.  [20]  incorporates  a dynamic  compensator  into  the  switching  function  (manifold) 
in  order  to  simplify  the  equivalent  control  law.  Then,  a control  law  that  is  a function 
of  the  switching  manifold  and  the  system  state  is  utilized  to  stabilize  a system  that 
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features  a constant  delay  in  both  a state  and  the  input.  The  work  by  Xia  et  al.  [55] 
considers  the  derivation  of  delay-independent  as  well  as  delay-dependent  stability 
conditions  for  a class  of  linear  systems  with  simultaneous  delay  in  the  state  and  the 
input.  An  integral  switching  function  with  a compensator  is  utilized  to  obtain  a 
simple  equivalent  control  law.  The  stability  conditions  are  given  in  terms  of  LMIs. 

This  chapter  introduces  a new  approach  to  address  the  problem  of  stabilizing 
a linear  system  featuring  both  state  and  input  delay.  First,  a state  transformation 
is  used  to  map  the  original  system  into  an  input-delay  free  form  where  only  state 
delays  are  present.  Then,  a new  state  defined  as  the  difference  between  the  original 
state  and  the  transformed  state  is  incorporated  into  the  transformed  system  equation. 
Introducing  an  integral  switching  function  in  terms  of  the  transformed  states  allows 
the  derivation  of  a simple  state  feedback  equivalent  control  law.  This  control  action 
along  with  a proposed  discontinuous  control  law  are  shown  to  drive  the  states  to 
the  switching  manifold  in  finite  time.  Finally,  using  a bound  on  the  new  state,  and 
utilizing  Lyapunov  techniques,  the  development  delivers  sufficient  stability  conditions 
that  ensure  the  asymptotic  stability  of  the  closed-loop  systems  in  terms  of  a constant 
LMI  that  depends  on  the  delay  of  the  system. 

The  chapter  is  structured  as  follows.  In  Section  3.2  the  problem  is  formulated 
along  with  the  transformation  that  eliminates  the  input  delay.  The  design  of  the 
control  law  which  consists  of  an  equivalent  control  action  and  a discontinuous  control 
action  is  discussed  in  Section  3.3.  In  section  3.4  the  control  law  is  shown  to  drive  the 
system  states  to  the  sliding  surface  in  finite  time.  Derivation  of  sufficient  conditions 
for  the  asymptotic  stability  of  the  transformed  and  original  systems  is  given  in  Section 
3.5.  A bound  on  the  time  delay  tolerable  by  the  systems  is  also  given.  The  chapter 
concludes  with  an  illustrative  example  that  verifies  the  results  in  Section  3.6,  and 
presents  a summary  in  Section  3.7. 
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3.2  Problem  Formulation 


The  time-delay  system  considered  is  of  the  form 


x{t)  = Ax{t)  + Adx{t  — h)  + Bu{t)  -f  Bdu{t  — h) 
x{t)  = $(r),  r€[— h,  0] 
u{t)  = ^'(r),  tG[— /i,  0] 


(3.1) 


where,  x{t)  G 5R"  is  the  state,  u{t)  G 3?"*  is  the  control  input,  and  A,Ad,B,  and 
Bd  are  matrices  of  appropriate  dimensions.  The  system  delay  h is  considered  to  be 
constant,  $(r)  is  an  initial-state  function,  and  ^(r)  is  an  initial-input  function.  The 
notation  | • | is  used  to  indicate,  dependent  on  the  scalar  or  vector  nature  of  the 
argument,  an  absolute  value  of  a scalar  quantity  or  a vector  norm,  and  |1  • ||  is  used 
to  indicate  an  induced  matrix  norm. 

System  Transformation.  The  following  state  transformation  is  introduced  as 
suggested  in  [1]  to  map  (3.1)  into  an  input-delay  free  system: 


Differentiating  equation  (3.2)  to  obtain 

z{t)  = x{t)  + A f e^^^-^-^^Bdu{r)dr  + e-^^Bdu{t)  - Bdu{t  - h) 

Jt-h 

and  then  substituting  for  x{t)  from  (3.1)  gives 

z{t)  = Az{t)  + Adx{t  — h)  + Bu{t) 

where  B = B+e~"^^Bd-  In  this  work  it  is  assumed  that  the  pair  {A,  B)  is  controllable. 
Let’s  define  the  new  state 


z{t)  — x{t)  + ( ^ '^^Bdu{T)dr 


(3.2) 


J t—h 


v{t)  :=  f ^ '^^Bdu{T)dr 


J t—h 


(3.3) 
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Then,  the  transformed  system  becomes 

z{t)  = Az{t)  + Adz{t  -h)  + Bu{t)  + Adv{t  - h)  (3.4) 

where  Ad  = -Ad-  Note  that  from  (3.2),  v{t)  = z{t)  - x{t)  is  interpreted  as  the 
difference  between  the  original  system  x{t)  and  the  transformed  system  z{t).  A 
feedback  matrix  F is  introduced  such  that  A = A- BF  is  Hurwitz  [47].  Treating  the 
last  term  in  (3.4)  as  an  internal  disturbance  and  defining  f{t,v{t-h)):=  Adv{t-h), 
the  system  equation  (3.4)  can  be  rewritten  as 

z{t)  = {A  + BF)z{t)  + Adz{t  -h)  + Bu{t)  + f{t,  v{t  - h))  (3.5) 

which  is  free  from  input  delay. 

3.3  Switching  Function  and  Control  Law  Design 
The  first  step  in  the  design  of  a sliding  mode  controller  is  to  define  a switching 
function  (manifold)  along  which  the  system  possesses  desired  properties,  such  as 
stability.  Various  structures  of  the  switching  functions  have  been  used  in  the  SMC 
literature.  The  most  common  designs  used,  however,  are  the  basic  form  s{t)  = Cx{t), 
the  integral  form  [47],  and  the  dynamic  or  compensated  form  [55].  The  basic  form 
is  best  suited  for  systems  having  the  general  structure  x{t)  = Ax{t)  + Bu{t).  The 
Integral  form,  adopted  in  our  work,  has  the  advantage  of  cancelling  the  delay  terms, 
which  allows  obtaining  a simple  state-feedback  equivalent  control  law  as  is  shown 
later.  The  dynamic  form  is  preferred  when  known  disturbances  and/or  delays  exist 
in  the  system,  as  it  helps  to  cancel  these  terms  hence  yielding  a simple  equivalent- 
control  law. 

The  sliding  surface  is  defined  by  a scalar  switching  function  s{t)  G 5R  of  the 


s{t)  = Cz{t) 


-1-  CAdz{r  — h)]  dr 


integral  form 


(3.6) 
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where  C is  a design  matrix  chosen  such  that  CB  is  nonsingular.  The  structure  of 
the  control  law  is  given  by 

u{t)  = Ue{t)  + Ud{t)  (3.7) 

where  Ue{t)  is  the  equivalent  part  and  Ud{t)  is  the  discontinuous  part  of  the  control 
law.  The  equivalent  control  is  obtained  by  setting  to  zero  the  derivative  of  equation 
(3.6)  with  respect  to  time,  and  then  solving  for  u{t)  to  yield 

s{t)  = Cz{t)  — CAz{t)  — CAdz{t  — h)  = 0 

Following  the  standard  approach  in  SMC,  the  state  derivative  z{t)  in  the  above 
equation  is  taken  from  (3.5)  after  ignoring  the  disturbance  term  — h)).  This 

gives  the  identity 

C{A  + BF)z{t)  + CAdz{t  -h)  + CBu{t)  - CAz{t)  - CAdz{t  - h)  = 0 
which  reduces  to 

CBFz(t)  + CBu(t)  = 0 

The  solution  to  the  above  identity  is  u(t)  = Ue(t);  hence,  after  recognizing  that  CB 
is  invertible  it  is  possible  to  conclude  that  the  equivalent  control  law  sought  is 

Ue(t)  = -Fz(t)  (3.8) 

The  discontinuous  control  law  proposed  is 

UdB)  = -(CB)~^[ks(t)  + p(t)  sgn(s(t))]  (3.9) 

where 

pit)  = \\C\\  ||i,||  \v{t-h)\  + C (3.10) 

and  where  A:  > 0 and  C > 0 are  design  parameters,  and  v{t—h)  = z{t—h)—x{t—h).  It 
must  be  noted  that  the  discontinuous  part  is  what  matters;  however,  the  linear  term 
in  equation  (3.9)  {i.e.,  ks{t))  helps  smooth  out  the  trajectories.  Various  structures 
of  the  discontinuous  control  law  can  be  found  in  Hung  et  al.  [26]. 
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3.4  Existence  of  a Sliding  Mode 

By  ” Existence  of  a sliding  mode”  we  mean  that  the  system  trajectories  must  be 
forced  to  reach  the  sliding  surface  in  finite  time  and  stay  there  forever.  Defining  a 
Lyapunov  function  V{t)  = ^s'^{t)s{t),  then  in  order  to  assure  reaching  the  manifold 
in  finite  time  it  suffice  to  show  that  V{t)  <0.  The  following  Theorem  provides  the 
proof. 


Theorem  1 The  time-delay  system  (3.5)  with  control  law  (3.7)-(3.10)  reaches  the 
sliding  manifold  within  a finite  time  tg,  where 


ts 


k |g(0)k 
C ^ 


(3.11) 


Proof:  Select  V{t)  - as  a candidate  scalar  Lyapunov  function.  Then, 

V{t)  = s{t)s{t)  = s{t){Cz{t)  - CAz{t)  - CAdz{t  - h)) 

= s{t){C{A  + BF)z{t)  + CAaz{t  - h) 

+ CB  [-Fz{t)  - {CB)~^{ks{t)  + p{t)  sgn{s))  ] 

+ Cf{t,  v{t  - h))  - CAz{t)  - CAdzft  - h)  } 

= s{t){-ks{t)  - p{t)  sgn{s))  + Cf{t,  v{t  - h)) 

= -k\s{t)\^  - p{t)\s{t)\  + Cf{t,v{t  - h))s{t) 

Now,  since  \Cf{t,v{t  - h))\  < HCH  ||id||  - ^)1,  after  invoking  (3.10)  it  follows 

that 

V{t)<-k\s{t)\^-C\s{t)\  (3.12) 


Therefore,  V{t)  < 0 for  all  /c  > 0 and  C > 0,  and  it  can  be  concluded  that 
the  system  trajectories  attain  sliding  mode  in  finite  time.  An  estimate  for  the 
upper  bound  of  the  reaching  time  tg  can  be  obtained  by  integrating  the  differen- 
tial equation  V{t)  = -k\s{t)\^  - Cls(i)|,  where  \s{t)\  = y/2V(t),  under  the  ini- 
tial condition  V"(0)  = |s(0)^.  The  result  (3.11)  is  obtained  after  a simple  trans- 
formation of  variables  and  a simple  algebraic  manipulations.  The  derivation  of 
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the  reaching  time  for  the  cases  k = 0 and  k ^ 0 is  presented  in  Appendix  A. 

■ 

3.5  System  Stability 

After  demonstrating  reaching  the  sliding  manifold  in  finite  time,  it  remains  to 
show  that  once  the  system  trajectories  are  in  the  sliding  phase  the  system  is  asymp- 
totically stable.  At  sliding  mode  the  control  law  (3.7)  reduces  to  u{t)  — Then, 
from  (3.8)  it  follows  that  the  dynamic  system  (3.5)  is  given  by  the  expression 

z{t)  = Az{t)  + Adz{t  - h)  + f{t,  v{t  - h))  (3.13) 

The  developments  in  the  suite  make  use  of  the  inequality 

\v{t)\  < rj{h)  \z{t)\  (3.14) 

where 

T]{h)  - h max  l|e~^®|l  11^^11  \\F\\a  (3.15) 

o<e<h 

and  where  a > 1 is  a constant  derived  using  a Razumikhin-like  argument.  The 
purpose  of  this  constant  is  to  describe  the  evolution  of  |2:(t)l,  i.e.,  \z{9)\  < o.  \z{t)\, 
0 € [t-h,t].  The  bound  (3.14)  follows  from  applying  successive  bounding  operations 
to  the  right-hand  side  of  (3.3)  and  introducing  the  Razumikhin  parameter.  The 
complete  derivation  of  the  constant  bound  in  (3.15)  is  given  in  Appendix  B. 

We  are  ready  now  to  provide  the  sufficient  conditions  for  the  asymptotic  stability 
of  the  system  (3.13)  which  are  introduced  by  the  following  theorem. 

Theorem  2 The  time  delay  system  (3.5)  with  control  law  (3.7)-(3.10)  is  asymp- 
totically stable  at  sliding  mode  if  there  exist  positive- definite  matrices  P G 
R G and  Q G 3?"^"  such  that 


^min{F)  > ^max  (Q) 


(3.16) 
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and 

Xrmn{Q){^rmn{R)  ~ XmaxiQ))  > (1  + VW?  \\PA£ 
where  P and  R are  solutions  to  the  Lyapunov  equation 

PA  + A^P  = -R 


(3.17) 


(3.18) 


Proof:  Consider  a Lyapunov  functional  of  the  form 

V{t)  = z^{t)Pz{t) [ z'^{T)Qz{T)dT  (3.19) 

Jt-h 

The  time  derivative  of  V{t)  with  respect  to  time  is  given  by 

V{t)  = 2z'^(t)Pz{t)  + {z^{t)Qz{t)  - z^(t  - h)Qz{t  - h)) 

Substituting  the  expression  for  z{t)  given  in  (3.13)  yields 

V{t)  = 2z^{t)PAz{t)  + 2z^{t)PAdz{t  - h) 

+2z^{t)Pf{t,  v{t  - h))  + z'^{t)Qz{t) 

—z^{t  — h)Qz{t  — h)  (3.20) 

and  then  using  (3.18)  and  bounding  the  right-hand  side  of  (3.20)  yields 

V{t)  < -Xrmn{R)\z{t)\'^  + ^max  {Q)\z{t)\'^ 

-\-2\\PAdW\z{t)\  \z{t-h)\ 

+2\\PAd\\  \z{t)\  \v{t  - h)\  - Xmin{Q)\z{t  - /i)|^ 

(3.21) 


Invoking  the  bound  (3.14)  and  rearranging  terms,  inequality  (3.21)  can  be  written  in 

the  form  ^ 

z(t)'^  a b z{t) 

V{t)  < (3.22) 


1 

T 

a b 

1 

■40 

z{t  — hY 

c d 

z{t  — h) 
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where 


a b 

(1  + vW)  \\PA,\\ 

c d 

(1  + ,(/.))  IIP.4,11 

Xmin  (Q) 

(3.23) 


To  prove  asymptotic  stability,  the  Lyapunov  functional  must  satisfy  V{t)  <0  which 
implies  that  the  matrix  (3.23)  must  be  negative-definite.  This  is  ensured  if  and  only  if 
conditions  (3.16)  and  (3.17)  are  satisfied.  B 

Theorem  2 can  be  reformulated  to  show  explicitly  the  constraint  on  the  size  of 
the  delay  parameter  imposed  by  design  choices,  such  as  the  adopted  Lyapunov  ma- 
trices R and  Q.  This  is  given  in  the  following  corollary. 


Corollary  1 The  time-delay  system  (3.5)  with  control  law  (3.7)-(3.10)  is  asymptot- 
ically stable  in  sliding  mode  for  time-delay  values  satisfying 

1 


h max  ||e  '^^||  < (-  — 1)t7j3  ,,  ,,  rp,, 

o<e</i"  ^2  \\PAd\\  ^\\Bd\\\\F\\a 


(3.24) 


where  P and  R are  positive-definite  solutions  to  the  Lyapunov  equation  (3.18)  and 
satisfy  the  inequality 


■^min(F)  > max(Aniax(Q)?  2||T*yld| 


(3.25) 


Proof:  The  proof  consists  of  deriving  conditions  that  ensure  the  existence  of  a feasible 
solution  to  (3.16)  and  (3.17).  The  proof  also  recognizes  that  XmaxiQ)  represents  the 
maximum  eigenvalue  to  the  Lyapunov  functional  (3.19),  hence  Xmin{Q)  < XmaxiQ)- 
Using  the  latter  inequality  along  with  the  constraint  imposed  on  XminiQ)  by  (3.17), 
it  follows  that 

A solution  XminiQ)  fo  (3.26)  exists  only  if 

il  + gih))VAdr 

XminiF)  Xmax  (Q) 


^ XmaxiQ) 


(3.27) 
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which,  using  the  fact  that  (3.16)  requires  that  Xmin{R)  ~ XmaxiQ)  > 0>  is  equivalent 
to 

XmaxiQ?  - X^UR)>^max{Q)  + (1  + V{h)nPA£  < 0 (3.28) 

The  analysis  of  the  above  inequality  reduces  to  investigating  the  boundary  defined 
by  the  equality 

XmaxiQ)^  - Xmin{R)Xmax{Q)  + (1  + T]{h)f\\PA£  = 0 (3.29) 

which  can  be  readily  solved  to  yield 

A™.(0)  = 4(1 + >,(/.) )^||Pyl,|P 

Given  that  only  real  solutions  are  meaningful,  it  follows  that  the  discriminant  must 
be  nonnegative,  i.e., 

Xrmn{R)>2{l  + rjih))\\PA,\\  (3.30) 

The  presence  of  the  factor  (1  + r]{h))  > 1 implies  that  a feasible  solution  to  (3.30) 
exists  only  if 

\min{R)  > 2\\PAd\\  (3.31) 

Furthermore,  from  (3.30)  it  follows  that  the  set  of  feasible  solutions  is  given  by  the 
equivalent  inequality 

T)(h)  < - — 1 

which  establishes  condition  (3.24)  of  the  corollary  after  using  (3.15)  and  suitably  re- 
arranging the  factors  in  the  inequality.  Moreover,  since  Xmin{R)  must  simultaneously 
satisfy  condition  (3.16)  and  constraint  (3.31),  it  follows  that  it  must  satisfy  condition 
(3.25)  of  the  corollary.  * 

It  remains  to  show  the  asymptotic  stability  of  system  (3.1),  as  addressed  in  the 


following  theorem. 
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Theorem  3 The  time-delay  system  (3.1)  with  state  x{t)  is  asymptotically  stable  if 
the  transformed  system  (3.5)  reaches  the  sliding  manifold  and  is  asymptotically  stable 
on  the  manifold. 

Proof:  If  z{t)  reaches  the  sliding  surface,  then  the  control  law  reduces  to  u{t)  - 
-Fz{t),  and  (3.2)  can  be  rearranged  in  the  form 


Now  when  (3.5)  is  asymptotically  stable,  it  follows  that  z{t)  — > 0,  then  x(t)  — ^ 0 in 
(3.32),  hence  completing  the  proof.  ■ 

The  ensuing  discussion  presents  an  example  that  illustrates  the  results.  From 
Figure  (3.1),  the  chattering  phenomenon  in  the  control  action  is  obvious.  Chattering 
refers  to  high-frequency  finite-amplitude  signals.  It  is  mainly  due  to  the  discontinuous 
control  law.  Chattering  is  undesirable  because  it  can  excite  neglected  high-frequency 
components,  and  lead  to  premature  wear  of  the  actuators. 

Several  approaches  have  been  proposed  in  the  literature  to  alleviate  or  eliminate 
the  effect  of  chattering.  Slotine  [49]  [50]  proposed  the  use  of  a boundary  layer  such 
that  standard  SMC  is  used  outside  the  boundary,  and  an  approximated  version  of  it 
takes  affect  inside  the  boundary.  The  work  of  Bartolini  [5]  introduces  a new  scheme 
to  chattering  reduction.  The  system  order  is  increased  and  an  estimator  based  on 
the  augmented  plant  is  defined.  Also,  a suitable  manifold  is  defined  such  that  the 
derivative  of  the  control  law  is  discontinuous  on  this  manifold.  Finally,  this  control 
law  is  feed  through  an  integrator  placed  in  the  plant  to  yield  a continuous  control  law. 
In  [34],  the  actuator  dynamics  are  treated  as  un-modelled  dynamics,  and  thus  are  not 
part  of  the  control  law.  Instead,  the  pass  filter  characteristics  of  the  actuators  are 
utilized  to  smooth  out  the  chattering  introduced  by  the  discontinuous  control  action. 


(3.32) 
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3.6  Example 

Consider  the  time-delay  system  (3.1)  with  h = 0.8,  with  an  initial-input  function 
^(r)  = 0 for  r G [—h,0),  and  an  initial-state  function  $(r)  = [—1  2]^  for  r 6 
[-h,  0]  so  that  the  initial-state  vector  is  x(0)  = [-1  2]^,  and  the  following  system 
parameters: 


-1  0 

0.01 

-0.04 

2 

1 

A = 

0.2  0.3 

, Ad  — 

0.02 

0 

, B = 

1 

, Bd  = 

0 

The  control  design  considered  is  based  on  the  following  matrices  associated  with 
the  Lyapunov  equation  (3.18)  and  the  Lyapunov  functional  (3.19): 


6 0 

3 0 

4.5435  5.7537 

R = 

, Q = 

, P = 

0 6 

0 3 

5.7537  48.4578 

Using  the  feedback  matrix  F = [-0.0166  0.2827]  the  eigenvalues  of  A are  placed 
at  {—0.4,  —0.45}.  The  controller  parameters  are  A:  = 30  and  C = 3-  The  switching 
function’s  initial  value  is  s(0)  = 3,  and  its  design  matrix  is  chosen  as  C = [0.9  0.85]. 
Selecting  a = 2.3391,  calculating  the  norms  of  \\Bd\\  = 1,  ||F||  = 0.2832,  and  eval- 
uating max|le“^®||  = 2.2381  , then  equation  (3.15)  gives  rj{h)  = 1.1858.  It  is  now 

0<9<h 

straightforward  to  verify  that  conditions  (3.16)  and  (3.17)  of  Theorem  2 are  satisfied. 
First,  condition  (3.16)  is  met  given  that  Xmin{R)  = 6 is  greater  than  XmaxiQ)  — 3. 
Also,  condition  (3.17)  is  met  since  Xmin{Q){Xmin{R)  ~ XmaxiQ))  = 9 is  greater  than 
{l  + r]{h))‘^  ||PA(i|p  = 5.4770.  It  follows  then  that  Theorem  2 ensures  the  asymptotic 
stability  of  the  closed  loop  system.  The  conditions  of  Corollary  1 are  also  satis- 
fied since  this  corollary  is  equivalent  to  Theorem  2.  In  fact,  = 1.7905 

is  less  than  " 2.7208,  and  A™„(F)  = 6 is  greater  than 

max(A^ax(Q),2||PAd||)  = 3. 
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Figure  3.1  shows  the  results  of  a simulation  study.  Figure  3.1(a)  depicts  the 
asymptotic  stability  of  the  transformed  system  (3.5)  with  state  variable  z{t).  The 
state  trajectories  for  the  original  system  (3.1)  with  state  variable  x{t)  are  shown  in 
Figure  3.1(b).  Equation  (3.11)  yields  tg  = 0.1145,  a value  that  is  consistent  with 
the  time  at  which  s{t)  becomes  identically  zero  in  Figure  3.1(d),  given  that  at  that 
instant  z{t)  has  reached  the  sliding  manifold.  It  is  apparent  that  the  states  x{t) 
develop  asymptotic  behavior  after  a time  t ^ tg  + h = 0.9145,  which  is  a consequence 
of  the  fact  that  the  original  system  has  an  input  delay  whereas  the  transformed  system 
is  free  of  input  delay.  Figure  3.1(c)  shows  the  control  action  u{t)  rising  quickly  from 
its  initial  value,  and  reaching  the  value  of  zero  at  approximately  the  same  time  that 
z(t)  reaches  the  sliding  manifold.  Figure  3.1(c)  also  shows  that  the  control  scheme 
suffers  from  a chattering  effect,  as  is  to  be  expected  from  the  presence  of  the  signum 
function  in  the  discontinuous  control  law  (3.9). 

The  chattering  of  the  signal  can  be  alleviated  by  introducing  the  approximation 

Figure  3.2  shows  that  a value  of  e = 0.001  effectively  makes  the  chattering  disappear 
(see  Figure  3.2(c)),  while  the  state  trajectories  z{t),  x{t),  and  the  switching  function 
s{t)  remain  virtually  unchanged. 

Remark  1 The  negative  definiteness  of  the  constant  matrix  (3.23)  can  be  checked 
directly  through  the  Linear  Matrix  Inequality  (LMI)  toolbox  of  Matlab. 

3.7  Conclusions 

A sliding  mode  controller  has  been  designed  to  stabilize  a linear  system  with 
state  and  input  delay.  A key  step  is  the  use  of  a transformation  to  map  the  system 
input-delay  free.  This  transformation  is  also  used  to  define  a new  state  appears  as 
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a disturbance  in  the  transformed  system.  A SMC  is  then  designed  to  stabilize  the 
state-delay  system.  The  controller  is  shown  to  successfully  drive  the  system  states 
to  the  sliding  surface  in  finite  time.  Sufficient  stability  conditions  are  derived  using 
Lyapunov  techniques. 


Figure  3.1:  Trajectories  for  the  states  of  the  transformed  system  (a),  the  states  of 
the  original  system  (b),  the  control  law  (c),  and  the  switching  function  (d). 
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Figure  3.2:  Trajectories  for  the  states  of  the  transformed  system  (a),  the  states  of 
the  original  system  (b),  the  control  law  (c),  and  the  switching  function  with  approx- 
imation to  the  signum  function  (d). 


CHAPTER  4 

STATE  FEEDBACK  CONTROL  OF  TIME-DELAY  BILINEAR  SYSTEMS 

4.1  Introduction 

This  chapter  considers  the  stabilization  of  a class  of  state-delayed  bilinear  system 
with  constant  delay.  Much  work  has  been  done  to  derive  sufficient  conditions  for  the 
asymptotic  stability  of  the  closed-loop  bilinear  systems  via  a variety  of  controllers, 
including  state  feedback  [24],  quadratic  feedback  [11],  and  nonlinear  control  [10,  12]. 
Work  has  also  been  done  to  derive  stability  conditions  for  time-delay  bilinear  systems. 
The  combination,  however,  of  the  time-delay  and  the  nonlinearity  makes  the  design 
of  stabilizing  controllers  as  well  as  the  analysis  much  more  challenging.  Stability 
conditions  can  be  either  delay-independent  or  delay-dependent.  Delay-independent 
conditions  do  not  give  any  information  regarding  the  size  of  the  delay  tolerable  by 
the  system  and  therefore  they  are  generally  more  conservative.  On  the  other  hand, 
delay-dependent  conditions  provide  information  about  the  bound  of  the  delay,  which 
leads  to  less  conservatism  [18].  The  former  provide  no  information  about  the  delay 
tolerable  by  the  system.  Some  results  can  be  found  in  [13,  23,  24,  42]  . In  Chiang 
[13],  the  stability  analysis  of  a class  of  input-delay  bilinear  systems  is  considered.  The 
derivations  utilize  the  Razumikhin  parameter  in  conjunction  with  matrix  measure 
techniques  [52].  The  stabilization  of  a class  of  state-delay  bilinear  systems  with 
saturating  actuators  is  investigated  in  Niculescu  et  al.  [42].  The  work  by  Ho  et  al. 
[24]  utilizes  a memory  state-feedback  control  law  to  yield  global  stability  conditions 
for  a class  of  time-delay  bilinear  systems.  In  Guojun  [23],  the  stabilization  of  a class 
of  time-varying  bilinear  systems  with  output-feedback  is  studied.  Delay-dependent 
conditions  are  given  in  Liu  [37],  where  a memoryless  state-feedback  control  law  is 
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used  to  derive  stability  conditions  for  a time-delay  bilinear  system  with  saturating 
actuators. 

In  this  chapter,  the  stability  analysis  of  a class  of  state-delay  bilinear  systems 
via  state  feedback  control  is  investigated.  Lemmas  2 and  3 are  developed  to  facilitate 
the  proof  of  the  main  theorem.  The  analysis  utilizes  the  matrix  measure  [52],  and 
a technique  that  allow  expressing  the  stability  conditions  in  terms  of  a bound  on 
the  system  delay  and  an  initial  condition  region  of  attraction.  As  a result,  delay- 
dependent  stability  conditions  are  derived. 

The  chapter  is  organized  as  follows.  Section  4.2  presents  the  problem  along 
with  a controllability  assumptions.  Section  4.3  of  the  paper  introduces  preliminary 
results  which  will  be  utilized  in  the  proof  of  the  main  result.  The  main  result  is 
presented  in  Section  4.4,  followed  by  an  example  and  conclusions  in  Sections  4.5  and 
4.6,  respectively. 


4.2  Problem  Statement 

Consider  the  system 

x{t)  = Ax{t)  -I-  Adx{t  - h)  + Bu(t)  -t-  Nx{t)u{t) 

x{(t>)  = ^(0),0G  [-h,0]  (4.1) 

where  t G 5R  is  the  time  variable,  x{t)  G 3?"  is  the  state,  u{t)  G 3?  is  a scalar  input, 
A,Ad,N,  and  B are  matrices  of  appropriate  dimensions,  and  is  an  initial-state 
function.  The  non-negative  system  delay  h is  considered  to  be  constant.  In  our 
development,  the  following  conventions  are  used:  the  notation  | • | is  used  to  denote  a 
vector  p— norm,  and  the  notation  ||  • ||  is  used  to  denote  the  induced  matrix  p— norm. 
Finally,  p(-)  is  used  to  denote  the  matrix-measure  function  (see  Appendix  E for 
definition  and  useful  properties)  based  on  the  induced  matrix  p-norm  . Also,  the 
following  two  assumptions  are  adopted: 
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(yll)  the  pair  {A^  B)  is  controllable. 

(A2)  the  pair  {A  + Ad,  B)  is  controllable. 

Our  objective  is  to  find  a linear  state  feedback  control  of  the  form 

u{t)  = Fx{t)  (4.2) 

that  renders  the  closed-loop  system  asymptotically  stable.  We  also  aim  at  deriving 
delay-dependent  stability  conditions  which  guarantee  the  asymptotic  stability  of  the 
system  for  initial  conditions  lying  in  a specified  region,  the  region  of  attraction.  Under 
the  feedback  control  (4.2)  the  closed-loop  system  becomes 

x{t)  = Ax{t)  -\-  Adx(t  — h)  + Nx{t)Fx{t)  (4.3) 

where  A = A + BF  is  Hurwitz  stable. 

4.3  Preliminary  Results 

The  main  stability  results  are  given  in  Theorem  1 for  the  proof  of  which  use  is 
made  of  the  following  three  lemmas. 


Lemma  1 Consider  the  scalar  differential  equation 


y{t)  =ay(tf  +by{t) 


(4.4) 


where  a > 0,  and  b ^ 0.  Then  the  analytic  solution  is  given  by  [13] 

fee*'*  ?/(0) 


y{t)  = 


where  y(0)  is  the  initial  condition. 


b + a 2/(0)  (1  — e***) 


(4.5) 


Proof:  The  derivation  of  the  analytic  solution  (4.5)  is  given  in  Appendix  C. 
Furthermore,  a thorough  discussion  of  the  solution  behavior  is  provided  in  Section  4.7 
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along  with  a graphical  interpretation  of  the  solution.  Now,  the  finite  escape  time,  tj, 
for  which  y{tf)  oo  can  be  found  by  setting  the  denominator  of  (4.5)  to  zero.  In 
our  development,  the  focus  is  on  non-negative  initial  conditions,  y{0)  > 0 that  are 
not  equilibrium  points  of  (4.4)  {i.e.,  y{0)  ^ 0 and  y(0)  “„)•  If  0 < y(0)  < — ^ and 

b < 0,  then  there  is  no  finite  escape  time,  and  from  (4.5)  the  following  observations 
are  readily  verified: 

fij  y{t)  > 0 V t < oo 

(ii)  y{t  + T)  < y{t)  V t < T < oo 

(in)  lim  y(t)  = 0 ■ 

t^oo 

Remark  1 It  should  be  noted  that  Claims  (i)  and  (Hi)  can  be  verified  from  the  ana- 
lytical solution  (4-5).  The  proof  of  Claim  (i)  is  given  in  Appendix  D.  Claim  (ii)  can 
be  proved  by  verifying  that  (4-5)  implies  the  inequality  y{t  -\-T)  < y{t)  V t < T < oo. 
The  complete  proof  is  given  in  Appendix  D. 

Lemma  2 Consider  the  scalar  differential  equation 

z{t)  = az{t)^  bz{t)  -I-  Co  (4.6) 

where  a > 0, 6 < 0,  and  Cq  > 0.  Let  ki  and  k2  be  the  roots  of  az{ty  bz{t)  -I-  Cq  . If 
5^  > 4aco,  and  ki  < ^(0)  < k2,  then 

(i)  z{t)  > 0 V t < oo 

(ii)  z{t  -t-  r)  < z{t)  y t <T  < oo 

(iii)  z{t)  — > /ci  as  t oo 

Proof:  The  condition  5^  > 4aco  implies  that  ki  and  k2  are  real  distinct  roots. 
Introducing  the  state  transformation 


v(t)  = 2(<)  + r 


(4.7) 
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where  r is  a real  constant  to  be  determined,  and  combining  the  time  derivative  of 

(4.7)  with  equation  (4.6)  yields 

y{t)  — a y{t)^  + {b  — 2ar)y{t)  + ar^  — br  + Cq  (4.8) 

The  quadratic  form  defined  by  the  last  three  terms  on  the  right-hand  side  of  equation 

(4.8)  is  set  to  zero  by  the  values 

b T Vb^  ~ 4aco  n\ 

ri,r2  = — (4.9) 

where  r\  < V2-  One  can  easily  verify  that  ki  = —r2  > 0 and  k2  — —r\  > 0.  The 
objective  is  to  select  a value  of  r that  renders  the  coefficient  b — 2ar  < 0 in  equation 

(4.8) .  This  is  realized  if  and  only  if  ^ < r,  which  as  a consequence  of  equation  (4.9) 
is  satisfied  by  selecting  r — r2,  where  T2  is  the  largest  root.  Therefore,  substituting 
r = T2  into  equation  (4.8)  yields 

y{t)  = a y(t)^  + {b  - 2ar2)y{t)  (4.10) 

Now,  equation  (4.10)  has  the  same  form  as  equation  (4.4);  hence.  Lemma  1 can  be 
applied  to  conclude  that  y{t)  ->  0 as  t — >■  oo,  provided  that  0 < y{0)  < — • From 

the  transformation  equation  (4.7),  it  follows  that  z{t)  —r2  = fci  as  t — > oo,  pro- 
vided that  the  initial  condition  satisfies  z(0)  < — — = k2-  This  proves  claim  (in)  of 
the  lemma.  Claims  (i)  and  (ii)  can  be  readily  verified  from  Lemma  1 and  the  transfor- 
mation (4.7).  First,  since  y{t)  > 0 then  z{t)  > —T2  > 0.  Second,  since  y{t)  is  strictly 
monotonically  decreasing  it  follows  that  its  shifted  version  z{t)  is  also  strictly  mono- 
tonically  decreasing.  ■ 

Remark  2 The  arrows  on  the  z-axis  of  Figure  4-io,  show  that  the  solution  converges 
to  the  smaller  equilibrium  point  ki  whenever  the  initial  condition  satisfies  z(0)  < k2- 
Figure  f.lb  shows  the  conceptual  state  trajectories  for  three  different  initial  condi- 


tions. 
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Figure  4.1:  Graphical  interpretation  of  the  differential  equation  of  Lemma  2:  (a) 
derivative  graph,  (b)  solution  curves 


Lemma  3 Consider  the  scalar  system 


z{t)  = az{t)^  + bz{t)  + hci  sup  z{9)  + hc2  sup  z{9)^  (4-H) 

t-2h<e<t  t-2h<e<t 


where  h > 0,  a > 0,  Ci  > 0,  C2  > 0, 5 < 0,  and  where  the  initial-state  function 


= z{(j))  >0,  (f)e  [-2/1, 0] 


(4.12) 


satisfies  sup  z{9)  = ^(0).  If 

-2h<0<O 


—b  — — 4aco 

2a 


< z(0)  < min{- 


b + hci 
a + hc2  ’ 


—b  + y/b"^  — 4aco 
2a 


} 


(4.13) 


and 

/i  < min{-^, — (4-14) 
4aco  Cl 

where  Cq  = h{ciz{0)  + C2z{0)'^),  and  Cq  = then 


(i)  z{t)  > 0 V / < oo 

(ii)  z{t  + T)  < z{t)  y t < T < oo 

(iii)  lim  z{t)  — 0 

t->00 

Proof:  Let  r = 2h.  Now,  consider  an  initial  condition  that  satisfies  (4.13),  hence 
2(0)  > 0,  and  assume  that  at  some  time  t2  < oo  the  state  function  satisfies  2(/2)  < 0, 


51 


and  that  the  state  changes  sign  for  the  first  time  at  an  instant  ti  < ^2  such  that 
h — t\  < T.  This  scenario  can  hold  only  if  z{t)  < 0 at  some  time  t G (^1,^2]-  The 
proof  consists  of  showing  that  this  is  a contradiction.  From  equation  (4.11)  it  follows 
that  for  all  t G (^1,^2]  the  state  derivative  z{t)  is  strictly  positive  because  all  the  term 
in  the  right  hand  side  are  positive.  This  contradicts  the  hypothesis  and  hence  proves 
claim  (i).  For  claims  (ii)  and  (in),  the  proof  is  conducted  in  four  steps.  First,  the 
system  is  shown  to  be  strictly  monotonically  decreasing  in  the  interval  [0,r).  Next, 
the  system  is  shown  to  be  strictly  monotonically  decreasing  in  the  interval  [r,  2r). 
In  a third  step,  it  is  shown  that  the  strict  monotonic  decrease  is  preserved  in  all 
subsequent  intervals  of  length  r.  Finally,  in  step  four  it  is  shown  that  a,st  ^ 00  both 
z{t)  — >•  0 and  z{t)  — >•  0. 

As  a preliminary  observation,  note  that  at  t = 0 equation  (4.11)  can  be  written 
as 

i(0)  = az{^Y  + ^-2^(0)  + hciz{0)  + hc2z{0Y 

= {a  + hc2)z{0)^  + {b  + hci)z{0)  (4-15) 

The  states  that  set  i(0)  = 0 can  be  found  by  finding  the  roots  of  the  expression 

0 = (a  + hc2)z(0)  [ 2^(0)]  (4-16) 

a + nc2 

Hence,  the  initial  states  z(0)  = Zi  and  2:(0)  = Z2,  where  = 0 and  Z2  = — a+hc2  P’-®" 
duce  zero  derivatives  at  the  initial  time  t = 0.  The  focus  now  turns  to  determining 
conditions  that  ensure  that  z{t)  is  a decreasing  function  of  time.  This  requires  that 
the  condition  i(0)  < 0 hold  at  all  finite  time. 

The  first  step  of  the  proof  considers  the  interval  t G [Or).  Depending  on  the  param- 
eter b,  there  are  two  scenarios  of  relevance.  If  6 + hci  > 0,  then  since  ^(0)  > 0,  it 
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follows  from  (4.15)  that  i(0)  > 0 and  the  solution  is  initially  increasing.  Obviously, 
this  case  is  not  desired.  However,  if 

b + hci<Q  (4.17) 


then  from  (4.15)  it  can  be  concluded  that  if 

^(0)< 


(4.18) 


a + hc2 

then  i(0)  < 0 and  it  follows  that  the  solution  is  initially  decreasing.  Therefore,  in 
the  interval  [0,r)  equation  (4.11)  can  be  written  in  a form  similar  to  that  given  in 
Lemma  2,  namely, 

Zo{t)  = azo{ty  + bzo{t)  + Co  (4-19) 


where  Cq  = hciz{0)  + hc2z{0)^  and  where  zo{t)  is  the  solution  of  (4.11)  for  t € [0, r). 
In  order  to  ensure  that  the  equilibrium  points  of  (4.19)  are  real  and  distinct  the 
discriminant  must  satisfy  5^  - 4aco  > 0,  which  implies  that  the  system  delay  must 
satisfy  h < hi  :=  Furthermore,  since  b + hci  < 0 is  the  desired  condition, 

it  follows  that  h must  also  satisfy  h < /12  :=  This  suggests  that  the  system 
delay  must  satisfy  h < min{hi,/i2}  which  yields  inequality  (4.14)  of  the  lemma. 
Let  ki  = ~'’~^2'a~^~  k2  = respectively  denote  the  smallest  and 

largest  roots  of  the  right  hand  side  of  (4.19).  Since  b < 0,  Lemma  2 can  be  applied 
to  (4.19)  to  conclude  that  the  solution  Zo{t)  is  strictly  monotonically  decreasing, 
and  that  zo{t)  ki  a,s  t ^ 00  provided  that  2:(0)  belongs  to  the  region  of  attraction 
ki  < 2;(0)  < /c2-  In  order  to  also  satisfy  the  constraint  (4.18),  the  region  of  converence 
is  redefined  as 

ki  < 2(0)  < min{-^~^^^\fc2}  (4.20) 

a + hc2 

which  is  equivalent  to  inequality  (4.13)  in  the  lemma. 
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The  second  step  aims  to  show  that  z(t)  decreases  in  [r,  2r).  Since  it  has  been 
established  that  z(t)  decreases  in  the  period  [0,r),  it  follows  that  sup  z{9)  = 

t<9<2t 

z{t  - t)  and  sup  z{6)‘^  = z{t  - r)^  in  the  period  [r,  2r).  Hence,  for  t G [r, 2r) 

t<6<2t 

equation  (4.11)  can  be  written  as 

Zr{t)  = aZr{t)'^  + bZr{t)  + c{t)  (4-21) 


where 

c{t)  = hciz{t  -t)  + hc2z{t  - tY  (4.22) 


and  where  Zr{t)  is  the  solution  of  (4.11)  when  t G [r,  2r).  Note  that  as  a consequence 
of  the  results  of  the  first  step  of  the  proof,  c{t)  is  strictly  decreasing  in  [r,  2r)  which 
implies  that  the  roots  of  the  right  hand  side  of  (4.21),  namely. 


Ki{t) 


—b  — y/b'^  — 4ac{t) 
2a 


(4.23) 


and  

-b  + - 4ac(() 

“ 2a 

are  such  that  the  smaller  root  Ki{t)  is  decreasing  and  the  larger  roots  K2Y)  is 
increasing.  This,  in  turn,  implies  that  z{t)  is  strictly  decreasing. 

The  third  step  involves  extending  the  results  of  the  second  step  to  the  subsequent 
intervals,  [2r,  3r),  [3r,  4r),  . . . , etc.  Let  z„^(t)  represent  the  solution  to  system  (4.11) 
in  the  interval  [nr,  (n  + l)r)  where  n > 2 is  an  integer.  When  t G [nr,  (n  + l)r) 
system  (4.11)  can  be  written  in  the  equivalent  form 


^miY)  — T bZji-j-{Y)  T c(t) 


where  c{t)  is  given  by  (4.22).  Note  that  when  n = 2,  the  function  c{t)  is  strictly 
monotonically  decreasing  in  t G [r,2r).  Repeating  the  argument  invoked  in  the 
second  part  of  the  proof,  namely  that  the  monotonicity  of  c{t)  in  the  interval  ensures 
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that  Ki{t)  given  in  (4.23)  is  strictly  monotonically  decreasing  in  that  interval,  leads 
to  the  conclusion  that  Znrit)  is  strictly  monotonically  decreasing  in  [nr,  (n  + l)r) 
when  n = 2.  The  proof  is  completed  by  induction  for  n = 3, 4,  . . . , etc.  Hence,  Znrit) 
is  strictly  decreasing  in  any  interval  of  length  r.  This  implies  that  z{t)  is  strictly 
decreasing  which  proves  claim  (ii). 

Step  four  of  the  proof  is  based  on  recognizing  that  from  claim  (i)  z{t)  is  bounded 
from  below,  and  using  the  fact  that  z{t)  is  strictly  monotonically  decreasing,  it  follows 
that  z{t)  -)•  L as  t ^ oo,  where  L < oo  is  a limit,  and  that  z{t)  ->  0.  Taking  the 
limit  as  t — )■  oo  on  each  side  of  equation  (4.21)  yields 

0 = (o  + hc2)L^  + (5  + hci)L  (4.24) 

Solving  for  the  two  limits  of  (4.24)  yields  Li  = 0 and  Now,  given  that 

the  initial  condition  (4.13)  implies  that  2:(0)  < — ^+^^2  “ as  t — > oo,  the 

decreasing  state  z{t)  must  reach  the  lower  limit  Li  = 0.  Thus,  z{t)  = 0.  This 

proves  claim  (in).  ® 


4.4  Main  Result 

Now,  utilizing  the  developments  in  Lemma  3,  we  are  ready  to  present  the  main 
result  which  concerns  the  asymptotic  stability  of  the  original  system  (4.1).  The 
approach  is  to  bound  the  norm  of  the  solution  of  (4.1)  (Le.,  ||a;(f)||)  by  a scalar 
function  z{t)  that  is  asymptotically  stable.  Thus,  when  z{t)  0 then  ||a;(t)||  — > 0. 
An  argument  based  on  the  comparison  theorem  [41]  (see  Appendix  F)  is  utilized. 


Theorem  1 The  time-delay  bilinear  system  (4-1)  under  assumptions  (Al)  and  (A2) 
and  the  state  feedback  control  law  (4.2)  is  asymptotically  stable  if 


0 < |x(0)|  < min{ 


b + hci 
a + hc2  ’ 


—b  + — 4aco 


} 


2a 


(4.25) 
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(4.26) 

4aco  Cl 

where  a = \\N\\  ||F||,  b = Ci  = ||^d^l|  + ||^dAdl|,  C2  = \\AdN\\  ||F|1,  Cq  = 

/ici|x(0)|  + hc2\x{0)\'^,  Co  = f,  and  A = A + Ad- 


Proof:  Since  equation  (4.3)  is  continuously  differentiable,  then 

ro 

x(t)  — x{t  — h)=  x(t  + 9)  dO  (4.27) 

J-h 

Substituting  for  x{t)  from  equation  (4.3)  and  rearranging  terms  yields  the  expression 
x{t  - h)  - x{t)  - { Ax{t  + 6>)  + Adx{t  + 9 - h)+  Nx{t  + 9)Fx{t  + 9)  }d9 

J-h 

which  can  be  substituted  in  equation  (4.3)  to  get 

i{t)  = Ax{t)+Ad  [a:(t)-  / { Ax{t+9)+Adx{t+9-h)+Nx{t+9)Fx{t+9)  }d9  ]+Nx{t)Fx{t) 

J-h 


X 


or 


/•o 

x{t)  = Ax{t)+  / (-Ad){  Ax{t+9)+Adx{t+9-h)+Nx{t+9)Fx{t+9)  }d9  +Nx{t)Fx{t) 

J-h 

(4.28) 

where  A = A + Ad-  The  solution  to  equation  (4.28)  has  the  form 

x{t)  = e^^xoF  [ [ {{-Ad)Ax{s + 9)  - Alx{s  + 9 - h) 

Jo  J-h 

-AdNx{s  + 9)Fx{s  + 0)}  d9  + Nx{s)Fx{s)  ] ds  (4.29) 


where  Xq  = a;(0)  is  the  initial  condition  obtained  by  setting  x(0)  = 4'(0)  in  equation 
(4.1).  Utilizing  the  matrix  measure  property  ||e^*||p  < [14]^  and  taking  the 
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norm  of  both  sides  of  equation  (4.29)  gives 

|x(t)|  < ^ f { \\AdA\\\x{s  + e)\  + \\Al\\\x{s  + e - h)\ 

Jo  J-h 

+\\AaN\\  1|F||||x(s  + 0)|M^^+ Ill'll  ll^ll 


Now,  the  inner  integral  can  be  bound  using  the  supremum  of  its  arguments  and  the 
length  h of  the  integration  interval,  to  yield 

|a:(t)|  < [ h \\AdA\\  sup  \x{e)\  + h\\Al\\  sup  la;(0)| 

Jq  s—hK9<s  s—2h<6<s~h 

+ h||/l,Afl|  ||F||  sup  |x(«)p  + ||iV||  ||F||  |x(s)|^  ] ds 

s—h<9<s 


Letting  a = ||A^||  HFjl,  ci  = \\AdA\\  + \\AdAd\\,  and  C2  = \\AdN\\  |1F||,  and  adopting 
the  largest  interval  s — 2h  < 9 < s gives 


\x{t)\  < e 


p.{A)t 


l^ol  + [ 

Jo 


=M(4)b-s) 


{ hci  sup  |a;(^)| 


s-2h<e<s 

+hc2  sup  |a;(0)p  + a|a;(s)p}ds 

s-2h<9<s 


(4.30) 


Now,  let  a scalar  function  with  initial  condition  Zq  — |a;o|  satisfy 

z{t)  = Zo+  [ { hci  sup  z{9)  + hc2  sup  z{9)^  + az{sf}ds 

Jo  s—2h<9<s  s—2h<9<s 

SO  that  \x{t)\  < z{t).  Then, 


z(t)  — az{t)‘^  + ii{A)z{t)  + hci  sup  z{9)  + hc2  sup  z{9Y  (4-31) 

t-2h<9<t  t-2h<9<t 

The  system  equation  (4.31)  is  of  the  same  form  as  that  of  equation  (4.11)  with 
b = IJ.{A).  Two  cases  are  considered,  depending  on  the  norm  of  the  initial  condition. 

Case  1:  xq  = x^  where  ki  < |a;ol  < min{-^±|^,  k2}. 

Case  2\  xo  = Xq  where  0 < |a:o|  < fci. 

where  ki  = -b-y/b^-iaco  ^ -b+y/b^-4aco  solution  to  (4.3) 

be  denoted  as  x'{t)  such  that  inequality  (4.30)  applies  with  x{t)  = x'{t).  Invoking 
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Lemma  3 show  that  z{t)  ^ 0 as  t ->  oo.  Hence,  since  \x'{t)\  < z{t)  it  follows  that 
x'{t)  — > 0 which  implies  that  x{t)  is  asymptotically  stable.  For  case  2,  let  the  solution 
to  (4.3)  be  denoted  as  x"{t)  such  that  inequality  (4.30)  applies  with  x{t)  = x"{t). 
Since  |a:d|  < |a;[)|  then  from  (4.30)  it  follows  that  \x"{t)\  < |x'(t)|.  Finally,  since 
x'{t)  ->  0 then  x"{t)  0,  which  implies  that  the  system  is  asymptotically  stable. 

This  completes  the  proof.  ■ 


4.5  Example 

Consider  the  time-delay  bilinear  system  (4.1)  with  h = 0.5,  an  initial-state 
function  4'(0)  = [—1.5  1.5]^  for  (j)  G [— 2/i,  0]  so  that  the  initial-state  vector  is 

x(0)  = [—1.5  1.5]^,  and  the  following  parameters: 


0.5  -0.2 

1 

O 

CM 

O 

1 

1 

-0.02  0.06 

) Ad  — 

, B = 

, N = 

0.8  -2.1 

0.5  -0.1 

0 

0.01  -0.03 

Choosing  the  eigenvalues  of  A to  be  placed  at  {—1.6  ,—2.1}  yields  the  feedback 
matrix  F = [-2.1  0.2].  Let  us  investigate  the  alternative  p-norms  to  verify  which 
norm  satisfies  the  inequality  conditions  (4.25)  and  (4.26). 

For  the  1— norm,  neither  condition  is  satisfied.  First,  the  discriminant  D = 
b‘^  — Aaco  = — 1.8  and  therefore  the  roots  are  complex.  Second,  /i  = 0.5  > {0.36, 0.06}. 
For  the  2— norm,  both  conditions  are  satisfied.  Hence,  adopting  the  2— norm  gives  the 
following  values  for  the  parameters  in  Theorem  1:  a = 0.1492, /i2(A)  = b = —1.32, 
and  Co  = 1.4130.  The  attraction  region  (4.25)  is  given  by  0 < |rc(0)|  = 2.12  < 
min{3.8,  7.6}.  Also,  the  bound  on  the  delay  is  given  by  h = 0.5  < min{1.18, 1.03}. 
Finally,  for  the  oo— norm,  both  conditions  are  satisfied.  The  values  of  the  parameters 
in  Theorem  1 are  given  by:  a — 0.168,  )Uoo(A)  = b = —1.8,  and  Cq  = 1.0415.  The 
delay-bound  condition  equation  (4.26)  is  satisfied  since  /i  = 0.5  < min{l. 44, 2.3148}. 
Furthermore,  the  domain  of  attraction  (4.25)  is  given  by  0 < |a;(0)|  = 1.5  < 

min{5.4855, 10.1}. 
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Figures  4.2  and  4.3  illustrate  the  results  of  a simulation  study.  Figure  4.2  depicts 
the  time  evolution  of  the  state  trajectories  which  assume  an  asymptotic  behavior. 
Figure  4.3  shows  the  norm  of  x{t)  converging  to  the  origin. 

Remark  3 The  sufficient  conditions  given  in  (4-25)  and  (4.26)  can  vary  depending 
on  the  norm  and  matrix  measure  chosen.  While  stability  can  be  concluded  for  a 
certain  norm,  it  may  not  be  so  for  other  norms.  The  stability  conditions  can,  however, 
be  tightened  and  therefore  reducing  conservatism  by  selecting  other  norms  and  matrix 
measures.  One  choice  could  be  the  following: 

II^IU  = max'^—  \aij\ 

3 

and 

\ ^ ^3  I n 

p,„,[A)  = max{au  + 2_^  — 10^1 1 
* ... 

4.6  Conclusions 

Delay-dependent  stability  conditions  are  derived  for  a class  of  time-delay  bilinear 
systems  utilizing  the  comparison  theorem  and  matrix  measure  techniques.  The  suf- 
ficient stability  conditions  provide  a bound  for  the  tolerable  system  delay  as  well  as 
the  domain  of  attraction  for  which  asymptotic  stability  is  guaranteed.  These  results 
are  however  conservative  mainly  due  to  applying  a supremization  over  a large  time 
interval  which  is  twice  the  size  of  the  delay. 
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Figure  4.2;  Plot  of  the  trajectories  of  the  system  states. 


Figure  4.3:  Plot  of  the  trajectories  of  the  system  states. 
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4.7  Further  Analysis  of  the  System  in  Lemma  1 
Further  analysis  of  the  system  equation  (4.4),  and  a discussion  of  the  behavior 
of  its  solution  (4.5)  is  presented  in  this  section. 

Theorem  2 Given  the  system  (4-4)  ^he  analytic  solution  (4-5),  then  the  follow- 
ing observations  are  readily  verified: 

Case  1:  b > 0 


OO 

if 

2/(0)  > 0 

(4.32a) 

lim  y{t)  = < 

t—^oo 

0 

if 

2/(0)  = 0 

(4.32b) 

_b 

a 

V 

if 

otherwise 

(4.32c) 

0 

if 

2/(0)  < 

(4.33a) 

lim  y{t)  = < 

t-¥ca 

_b 

a 

if 

11 

1 

(4.33b) 

OO 

if 

otherwise 

(4.33c) 

Proof:  For  Case  1,  it  is  noted  from  (4.5)  that  the  solution  becomes  unbounded 
at  finite  time  when  the  denominator  is  zero.  The  finite  escape  time  (FET)  can  be 
calculated  as  a function  of  the  parameters  a,b,  and  y(0),  by  setting  the  denominator 
to  zero  and  solving  for  t.  The  result  is 

Ifet  = T (4-34) 

b ay{0) 

Now,  to  prove  the  first  branch  of  Case  1,  namely  (4.32a),  the  denominator  of  (4.5)  is 
set  to  zero,  i.e., 


5 + ay(0)(l-e“)  = 0 
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or 

1 + ^ y{0)  (1  - = 0 (4.35) 

Now,  since  as  t — )•  oo  the  quantity  1 — e***  takes  values  between  zero  and  — oo,  and 
since  f?/(0)  > 0,  then  the  equality  (4.35)  is  satisfied  for  some  t > 0.  Hence,  the 
denominator  becomes  zero  and  the  solution  blows  up  at  finite  time.  The  second 
branch  of  Case  1 (4.32b)  is  readily  verified  by  substituting  y{0)  = 0 in  (4.5).  Finally, 
substituting  for  y(0)  = into  (4.5)  gives  y{t)  = which  proves  the  branch 
(4.32c).  Graphs  (a)  and  (b)  of  Figure  4.4  show  the  derivative  graph  and  the  solution 
curves  of  (4.4),  respectively,  for  the  case  where  b > 0. 


Figure  4.4:  Graphical  interpretation  of  the  differential  equation(4.4)  for  the  case 
where  b > 0 : (a)  derivative  graph,  (b)  solution  curves 

For  Case  2,  the  first  branch  (4.33a)  indicates  that  the  system  is  asymptotically 
stable  when  the  initial  condition  satisfies  ^(0)  < — Therefore,  in  this  case  the  finite 
escape  time  is  avoided.  To  verify  this  claim,  consider  the  denominator  expression 
(4.35).  For  t > 0 the  quantity  1 - e***  takes  values  in  the  interval  (0  , 1).  The 
goal  is  to  show  that  for  ^(0)  < the  equality  (4.35)  is  never  satisfied.  First,  it  is 
trivial  to  verify  this  claim  for  y(0)  = 0.  Second,  for  0 < y{0)  < — Mt  follows  that, 
since  — 1 < | y(0)  < 0,  the  second  term  in  (4.35)  is  never  equal  to  1.  Finally,  for 
y(0)  < 0,  equality  (4.35)  is  again  never  satisfied  since  its  second  term  in  is  always 
positive.  This  proves  the  first  branch  of  Case  2.  The  second  branch  (4.33b)  is  readily 
verified  by  substituting  y{0)  ^ in  (4.4).  Finally,  branch  (4.33c)  indicates  that 
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for  y{0)  > —^  there  is  a finite  escape  time,  which  means  that  (4.35)  is  satisfied  for 
some  finite  values  of  t.  First,  note  that  -1  < f y(0).  Next,  since  the  values  of  the 
quantity  1 - e***  belongs  to  the  interval  (0  , 1),  then  for  some  time  t > 0 the  second 
term  in  (4.35)  can  equal  -1  such  that  (4.35)  is  satisfied.  This  completes  the  proof 
of  Case  2.  Figure  4.5  depicts  the  solution  behavior  for  the  case  where  b < 0.  It  is 
clear  that  when  y(0)  < the  solution  converges  asymptotically  to  the  origin,  and 
for  y(0)  > —^  the  system  becomes  unstable. 


Figure  4.5:  Graphical  interpretation  of  the  differential  equation(4.4)  for  the  case 
where  b < 0 : (a)  derivative  graph,  (b)  solution  curves 


CHAPTER  5 

FUTURE  WORK  AND  DISCUSSIONS 


Future  work  may  focus  on  three  problems.  First,  designing  a sliding  mode 
controller  for  a class  of  time-delay  bilinear  systems.  Second,  extending  the  concept 
of  the  Nyquist  robust  sensitivity  margin  (NRSM)  to  a class  of  uncertain  systems  with 
multilinear  uncertainty  structure.  Finally,  the  concept  of  the  NRSM  can  be  utilized 
to  obtain  a weighting  function  for  an  design  similar  to  work  by  Baowei  [?]. 

5.1  Problem  1 : Sliding  Mode  Control  for  a Delayed  Bilinear  System 

The  system  under  consideration  is  of  the  form 

x{t)  = Ax{t)  + Adx(t  -h)-\-  Bu{t)  + Nx{t)u{t) 

x{t)  = $(r),  r€[— h,  0]  (5.1) 

where  the  system  delay  h is  assumed  constant. 

The  objective  is  to  design  a sliding  mode  control  (SMC)  that  renders  the  sys- 
tem asymptotically  stable.  However,  when  nonlinearity  is  combined  with  time-delay, 
the  problem  becomes  more  challenging  from  the  control  design  viewpoint  as  well  as 
from  the  perspective  of  the  stability  analysis.  Therefore,  problems  may  arise  when 
designing  a controller  to  force  the  trajectories  into  the  sliding  manifold  and  to  ensure 
they  stay  there  for  all  subsequent  time. 

5.2  Problem  2 : Extending  the  NRSM 

Multilinear  Uncertainty.  An  interesting,  yet  challenging,  problem  is  to  ex- 
tend the  concept  of  the  Nyquist  robust  sensitivity  margin  to  include  linear  system 
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with  multilinear  uncertainty  structure,  i.e.,  systems  of  the  form 


(5.2) 


where  r and  q are  multilinear  vectors.  In  the  polynomial  case,  since  the  multilinear 
uncertainty  lacks  edge  results,  the  mapping  theorem  [7]  introduces  overbounding 
polynomials  which  of  course  can  be  conservative.  For  the  transfer  function  (5.2), 
there  are  two  situations.  First,  when  the  vectors  r and  q are  independent,  stability 
can  be  analyzed  using  the  mapping  theorem  by  considering  a polytopic  family  g{s) 
such  that  any  worst-case  margin  calculated  for  g{s)  is  considered  a guaranteed  margin 
for  g{s)  [7].  However,  when  the  vectors  r and  q depend  on  each  other,  then  there  are 
no  comparable  results  to  analyze  the  robust  stability  of  the  system. 

H -oo  Design.  Another  future  work  project  is  to  design  an  H-oo  controller  based 
on  a weight  function  derived  from  A)jv,s-  This  project  can  follow  the  work  of  Baowei 
[29,  30]  which  is  summarized  as  follow.  Given  the  system  in  Figure  5.1,  use  of  made 
of  the  M - A structure  given  in  Figure  5.2  where  it  is  known  from  the  small  gain 
theory  [57]  that  the  system  in  Figure  xx  is  internally  stable  for  any  A(s)  satisfying 


r 


u 


Figure  5.1:  The  negative  feedback  loop  of  the  uncertain  system  p{s)  with  a controller 
c(s). 
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The  transformation  of  the  system  in  Figure  5.1  into  the  M — A formulation  is 
given  in  Figure  5.3  where  it  follows  that  the  system  is  stable  for  any  5(s)  satisfying 


||5(s)|loo  < 


1 


(5.3) 


where  R{s)  = Next,  the  system  in  Figure  5.3  is  put  into  the  mixed  sensi- 

tivity framework  where  the  stability  conditions  can  be  expressed  by  the  inequality 


lllF2(s)i?(s)||oo  < 1 


(5.4) 


The  problem  now  is  to  choose  W2{s)  to  represent  the  effective  part  of  S{s).  Fi- 
nally, a weighting  scheme,  namely  the  effective  critical  perturbation  radius  (ECPR), 
is  designed  based  on  the  critical  direction  theory. 


Figure  5.3:  A system  with  parametric  uncertainty  in  the  standard  M — A loop. 


APPENDIX  A 

DERIVATION  FOR  THE  REACHING  TIME 


General  case 

Consider  a Lyapunov  function  V{t)  — The  standard  condition  that 

ensures  reaching  the  sliding  manifold  in  finite  time  is  given  by 


V{t)  = s{t)s{t)  < -?7|sWI 


(A.1) 


where  t]  > 0 [49]. 


Theorem  1 Under  the  inequality  given  in  (A.l),  the  time  at  which  the  sliding  man- 
ifold is  reached  is  given  by 

ts  - -|so|  (A.2) 

V 

Proof:  The  proof  utilizes  Figure  A.l.  From  the  figure  it  is  noted  that  the  initial 
value  of  the  sliding  function  is  s(0),  and  that  s{ts)  = 0.  The  equality  limit  of  (A.l) 
(f.e.,  s{t)s{t)  — — 77|s(t)|  ) can  be  written  as 


dt  = — 


77|s(t)| 


ds 


(A.3) 


where  s{t)  ^ 0.  The  reaching  time  tg  is  obtained  by  integrating  (A.3)  as  follows: 


which  yields  the  result 


I -J(s(^s)  - s(0))  : s(t)  > 0,  case{a) 

I -^(-s(ts) + s(0))  : s(t)  < 0,  case{b) 


(A.4) 
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Since  s{ts)  = 0,  equation  (A.4)  can  be  written  as 


ts  = 


^(0) 


or 

ts  = ^ls(0)| 


40 


s{t) 


Figure  A.l:  Plot  of  the  switching  function 


Case  (b) 

for  s(0)  >0  (a),  and  s(0)  <0  (b). 


New  condition 

Consider  the  Lyapunov  condition 

y(t)  = -p|sp  - d|s|  (A. 5) 

where  p > 0 and  d > 0.  Using  the  fact  that  \s{t)\  = vW,  equation  (A.5)  can  be 
rewritten  as 

V{t)  = -2pV{t)  - y/2dV{tf^  (A.6) 

Let  y{t)‘^  = U(t).  Then  equation  (A.6)  becomes 

V{t)  = 2y{t)y{t)  = -2py{tf  - V2dy{t)  (A.7) 


or 


y{t)  + py{t)  = -j^d 
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Multiplying  through  by  and  rearranging  terms  yields 


Integrating  both  sided  of  (A. 8)  gives 


y{t)eP^  = — e'’*  + c 

The  constant  c can  be  found  by  evaluating  (A. 9)  at  t = 0,  which  gives 

^ = 2^(0)  + Wp 

Substituting  for  the  constant  c,  equation  (A.9)  can  be  written  as 

y{t)  = e“'’*[y(0)  + 


\/2p  \/2p 


Now,  using  (A. 6)  equation  (A.  10)  can  be  written  as 

v/V(i)  = e-o'iym  + 

Since  at  sliding  mode  V(t)  — yJV{t)  = 0 because  s{t)  = s{t)  = 0,  it  then 
that 

0 = -pt,  + ln\yW)  + ^]-\n^^ 


or 


, 1 , 

*S  = - In  1 H 1 

p a 


Now,  since  V{t)  = |s(t)^s(t),  it  follows  that 


1, 

tg  — — ln[l  H 


s(0)^s(0) 
2 


(A.8) 


(A.9) 


(A.10) 


follows 


(A.11) 


d 


APPENDIX  B 

DERIVATION  OF  THE  BOUND  ON  v{t)  OF  CHAPTER  3 


This  appendix  presents  a derivation  of  the  bound  (3.14) 

\v{t)\  < v{h)\z{t)\ 


Consider  the  equation 

v{t)  = [ (B.l) 

J t—h 

Substituting  for  the  control  law  u{t)  = -Fz{t)  valid  at  sliding  mode  and  taking  the 
norm  of  both  sides  of  (B.l)  yields 


k(i)|  < [ ^ ^^11  \\Bd\\  \\F\\\z{T)\dr  (B.2) 

Jt-h 

Let  0 = -(t  - h - r).  Then  for  r = t,  0 = h,  and  for  r = t - h,  0 = 0.  The  inequality 
(B.2)  then  implies 


where  z{t,  t — h) 


k(i)|  < f l|e  \\Bd\\  \\F\\^maxJz{i;)\de 

JO  — 

< h maxlle^"^®  \\Bd\\  ||Fl|  z{t,t-h)  (B.3) 

= max  \z(ib)\.  Utilizing  the  Razumikhin  concept  [38],  it  follows 


that 


z{t,  t — h)  < a \z{t)\ 

where  a is  a Razumikhin  parameter.  Therefore,  inequality  (B.3)  can  be  written  as 


\v{t)\  < h max  ||e-^'l|  ||Rd||  |1F||  a\z{t)\  (B.4) 

OKd^h 

which  yields  the  bound  |u(t)l  < r]{h)\z{t)l  where  r]  represent  the  right-hand  side  of 
inequality  (B.4). 
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APPENDIX  C 

SOLUTION  FOR  A DIFFERENTIAL  EQUATION 

Here,  the  derivation  of  the  solution  (4.5)  of  the  system  equation  (4.4)  given  in 
Lemma  1 of  Chapter  4 is  presented. 

Theorem  1 Consider  the  scalar  differential  equation 


x(t)  — ax{tff  + bx{t) 


(C.l) 


then  the  analytic  solution  is  given  by 


x{t) 


fee***  x(0) 


b + a a;(0)  (1  — e^*) 
Proof:  Equation  (C.l)  can  be  written  as 

dx 


Integrating  both  sides  gives 


/ 


bx{t)  + ax{tff 


dx 


= dt 


— t p C\ 


(C.2) 


(C.3) 


x{t){ax{t)  + b) 

where  ci  is  the  integration  constant.  Working  out  the  integral  yields 

1 . r 1 . 

— Inf 7— r r]  — t Cl 

b ^ax{t)  + b 

which  after  multiplying  through  by  b,  and  taking  the  exponential  of  both  sides  gives 


= e''*C2 


ax{t)  + b 


(C.4) 
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where  C2  = At  i = 0,  the  constant  C2  = Substituting  back  into  (C.4), 

and  re-arranging  terms  yields 


x{t)[ax{0)  + b]  = e’’^x{Q)[ax{t)  + b] 


Finally,  solving  for  x{t)  gives  the  solution 


x{t) 


be’’*  a;(0) 


6 -f  a x(0)  (1  — e***) 


APPENDIX  D 

PROOF  OF  CLAIMS  (i)  & (ii)  OF  LEMMA  1 OF  CHAPTER  4 


For  convenience  we  rewrite  the  system  equation  of  Lemma  1.  Given  the  scalar 
differential  equation 

y{t)  = ay{tf  + hy{t)  (D.l) 


where  a > 0 and  b < 0.  Its  analytic  solution  is  given  in  Appendix  C and  can  be 
written  as 


y{t)  = 


e*'*  y{0) 


1 + f y(0)  (1  - eW) 


(D.2) 


Now,  for  the  initial  condition  bound  0 < y(0)  < we  need  to  prove  the  following 
two  claims: 

Claim  (i)  : y{t)  > 0 V t < oo. 

Claim  (ii)  : y{t  + T)  < y{t)  V t < T < oo. 

Claim  (i)  can  be  checked  by  verifying  that  both  the  numerator  and  the  denom- 
inator of  (D.2)  are  either  positive  or  negative.  Let’s  consider  the  numerator  first. 
Since  the  exponential  function  is  always  positive  (z.e.,  e***  > 0),  and  y(0)  > 0,  then 
the  numerator  is  positive.  Now,  for  the  denominator,  since  0 < y(0)  < -^  then  it 
can  be  readily  verified  that 

-1  < -MO)  < 0 (D.3) 

b 

Furthermore,  using  the  fact  that  0<1  — e*’^<lVt>0it  follows  that 


0 > ^ y(0)  (1  - > -1 


(D.4) 


Therefore,  the  denominator  of  (D.2)  is  easily  seen  to  be  positive.  Hence,  y{t)  > 
OV  t < oo.  This  completes  the  proof  of  claim  (i). 
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Claim  (ii)  implies  that  the  solution  is  strictly  monotonically  decreasing.  It  can 
be  checked  by  verifying  that  the  ratio  i?  < 1,  where 

y{t  + T) 


R 


2/W 


l + fy(0)(l-e''{‘+^)) 


obt 


which  after  simple  manipulations  reduces  to 

1 + |y(0)(l  — 


R 


(D.5) 


is  satisfied  for  all  t < T < oo.  Using  (D.2),  the  ratio  (D.5)  can  be  written  as 


(D.6) 


(D.7) 


Multiplying  the  numerator  and  the  denominator  of  (D.6)  through  gives 

gftr  ^ e'’^fy(0)  - e^^e'’*fy(0) 

1 + fy(0)  - e'’^e^*f?/(0) 

Ignoring  the  last  term  in  the  numerator  and  denominator  of  (D.7)  since  they  are  the 
same,  then  the  truncated  ration  R!  is  given  by 


„ e‘^(l  + Ivm 

i + !»(o)  ■ 


(D.8) 


for  all  r > 0.  This  shows  that  i?  < 1,  and  hence,  proves  claim  (ii)  of  Lemma  1. 


APPENDIX  E 
MATRIX  MEASURE 


The  following  definition  of  the  matrix  measure  and  its  properties  are  found  in 
Vidysagar  [52]. 


Definition  1 The  matrix  measure,  also  known  as  the  logarithmic  derivative,  of  an 
induced-matrix  norm  ||  • ||p  on  is  a function  fip  : — > TZ  defined  by 


Hp{A)  = lim 

e-yoo 


||J  + eA||p  - 1 
e 


(E.l) 


The  matrix  measure  of  A G corresponding  to  the  1,2,  and  oo  norms  is 

given,  respectively,  by 

//pi (A)  = max{ajj -\-y  \apj\} 

hp2{A)  = \jnax[{A*  A)\l 2 
/^poo(A)  — max^Opp  T ^ 

Some  useful  properties  of  matrix  measure  include  the  following: 

• /ip(A  B)  = fj-p(A)  + fJ^p{A). 

• —iXp{—A)  < Re  X < //p(A)  where  A is  an  eigenvalue  of  A. 

• //p(-)  is  a convex  function  on 
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APPENDIX  F 

THE  COMPARISON  THEOREM 


Let  a vector-valued  function  v{t,  s,z)  : J x J x > 5?"*,  J :=  [to,  oo]  has  the 
following  property.  For  any  fixed  t,  s, 

Zi  < Z2^  V{t,  S,  Zi)  < V{t,  S,  Z2) 

Let  z{t)  be  the  solution  to  the  inequality 

< -z(O)  + f v{t,s,z{s  — h))  ds 
Jto 

Then  the  maximal  solution  r{t)  of 

w{t)  < 2:(0)  + f v{t,  s,  w{s  — h))  ds 

Jto 

satisfies  z{t)  < r{t)  for  t > to- 
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